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ABSTRACT 


Conventional pitot-static airspeed measurement systems do not yield accurate 
measurements when aircraft speed is below 40 knots. Recent studies have demonstrated 
that neural network approaches for predicting airspeed are quite promising. In this thesis, 
a back-propagation neural network is used to predict the airspeed of UH-60A and OH-6 A 
helicopters in the low speed environment. The input data to the neural networks were 
obtained using the FLIGHTLAB flight simulator. The results obtained by flight 
simulation were validated by comparison to results of a previous study of the UH-60A 
helicopter based on actual flight data. The results of the work performed for this thesis 
show that at sea level the UH-60A low airspeed can be predicted with an accuracy of 
±0.71 knots and ±0.88 knots for out of ground effect and in ground effect conditions 
respectively. OH-6A analyses were performed at two pressure altitudes. At sea level the 
OH-6A airspeed can be predicted with an accuracy ± 0.75 knots when the aircraft is out 
of ground effect and ±0.88 knots when the helicopter is in ground effect. At a pressure 
altitude of 6000 feet OH-6A airspeed can be predicted with an accuracy of ± 0.64 knots 
for both flight conditions. 
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I. INTRODUCTION 


Since the 1970s there has been a signifieant eoneern to inerease the aeeuracy of 
airspeed measurement systems beeause eurrent measurement systems are ineffective at 
low airspeeds. The focus of this thesis is to develop a neural network (NN) model to 
estimate helieopter airspeed in the low speed environment using the NeuralWorks 
Professional Plus/ll software. A flight simulator provides the required data that is used as 
inputs to the neural network model. 

The thesis is organized into five chapters. Chapter 1 is the introduction chapter, 
whieh presents research results eurrently available in the area of neural network and 
aireraft speed measurement systems. Chapter 11 introduees the main eoneepts behind NN 
implementations and the basie strueture of the Baek-Propagation Neural Network 
(BPNN). Chapter 111 deseribes a speeifie implementation using the NeuralWare 
Professional Plus/ll software and presents the FLIGHTLAB simulation tool. The 
simulations performed as part of this thesis are diseussed in Chapter IV. Finally, Chapter 
V summarizes results obtained and discusses avenues for further researeh. 

A. BACKGROUND 

Sinee 1995, the Naval Postgraduate Sehool (NPS) has eondueted several researeh 
projeets using OH-6A “Cayuse” helieopters after reeeiving two sueh helieopters from the 
Massaehusetts Army National Guard. One of these projeets involves the development of 
a Vortex-Ring State Warning System (VRSWS). The vortex-ring state is a eondition of 
powered flight where the helieopter settles into its own downwash. The helieopter will 
inerease its rate of descent very rapidly as the lift is lost when entering the vortex-ring 
state, and any further applieation of eolleetive, a flight eontrol meehanism for helieopters, 
tend to reduce rotor effieieney. In this state the rotor experiences a very high vibration 
level and loss of eontrol. The eonsequenees of the vortex-ring state when the helieopter is 
elose to the ground might be extremely dangerous beeause loss of eontrol at low altitudes 
often results in aireraft erash. Development of VRSWS requires better airspeed 
measuring systems than these eurrently available with most avionies systems. 
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Last year, NFS signed a Cooperative Researeh and Development Agreement 
(CREDA) with Advanee Rotoreraft Teehnology, Ine. (ART) to reeiprocally develop 
advaneements in rotoreraft teehnology. Based on this agreement, NFS provided ART 
with equipment for a flight simulator, including the sticks and grips, seats and avionics 
suite of an OH-6A for integration into the control loader platform. ART provided their 
helicopter modeling suite with advanced visual rendering equipment to produce a fully 
functional stationary open platform flight simulator. In return for this equipment, ART 
provided the flight modeling software to NFS to develop rotoreraft models for future use 
which includes the OH-6 A and V-22 “Osprey” models. NFS was able to obtain flight 
parameters such as rotor RFM, engine torque, roll rate, etc., using simulation techniques 
through its cooperation with ART. LT Gregory Ouellette is currently conducting research 
about FLIGHTLAB, a flight simulation tool provided by ART. The required flight 
parameters for the NN model of this thesis were provided by LT. Ouellette, USN 
[Ref 19]. 


B. TECHNOLOGICAL PROBLEM 

The aerodynamic velocity in the plane of a helicopter’s rotor disc significantly 

affects the control characteristics of that helicopter. The pilot easily feels this effect 

during flight and especially in hover, low speed and transition conditions. Small changes 

in these flight regimes are rapidly followed by substantial changes in control sensitivities 

and trim positions. Therefore, the pilot’s job is very difficult when performing tasks 

during these specified conditions, and it is possible for these conditions to occur in any 

flight mission. Filots are capable of successfully handling the control helicopter with 

training but in IFR flight, where ground references are no longer accessible, the pilot’s 

job becomes increasingly more difficult. The need to extend military operations during 

poor weather and at night has been inevitable for military effectiveness throughout 

history. Therefore, many military missions have been performed during poor visibility 

and in a low speed environment. However, operations under these conditions endanger 

flight safety. The roles of search and rescue, submarine detection, mine countermeasures, 

front line supply, air-to-ground attack, and reconnaissance all demand a low speed, poor 

weather, day or night capability. Thus, several technological improvements have been 
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made to aircraft to conduct missions without endangering the lives of the aircrew. Today, 
military helicopters perform a wide variety of tasks in conditions ranging from hot and 
dry to cold and wet, windy and low visibility weather. Accurate low speed velocity 
sensing devices are essential because aircraft velocity and position information are what 
pilots need to perform the aforementioned tasks. 

Conventional methods to measure airspeed for aircrafts have been used for over 
60 years. Pitot-static systems are still commonly used since they offer simple, low cost 
and reliable enough solutions to measure airspeed. In this system, airspeed is derived by 
measuring the difference between the total pressure occurring at the forward facing pitot- 
probe and the static pressure measured at a static vent [Ref. 12]. Since helicopters fly at 
relatively lower speeds than aircrafts, and the cruise speed of a helicopter is less than 0.3 
Mach, helicopters are used to fly in incompressible subsonic flow conditions where 
Bernoulli principles pertain. In such flow regimes, airspeed is obtained simply by taking 
the square root of this pressure difference and multiplying the scale factor. 

The flight of a helicopter occurs in the two distinct regimes of hover/low speed 
flight up to 45 knots, including vertical maneuvering, and mid/high speed flight up to Vne 
-never exceed velocity- where Vne is the maximum airspeed permitted under any 
circumstances [Ref 17]. It is defined as a function of altitude, temperature and gross 
weight. For example, using flight manual of UH-60A, Vne is found to be 186 knots at -20° 
C, 4000 feet and when the gross weight of helicopter islSOOO lb. The low speed regime is 
very much unique to the helicopter as an operationally useful regime. No other flight 
vehicles are as flexible and efficient at maneuvering slowly close to the ground and at 
avoiding obstacles. Therefore, the low speed regime is a significant portion of helicopter 
flight time. In fact, the maneuvers done in this regime make the helicopter invaluable. 
Although low speed is very critical for helicopters during take-off, landing or hovering, 
measuring airspeed and wind direction is generally lacking in this regime [Ref. 3]. 
During low speed flight, the current airspeed measurements systems are inaccurate due to 
the rotor downwash and limitations of the pitot-static system. The conventional pitot- 
static sensor is ineffective at airspeeds below 40 knots and does not function at all during 
rearward flight [Ref. 8]. John Carter explains in Ref. 12 why it does not function 
properly. 
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The pressure dijference in pitot-static system becomes very small at low 
airspeeds. This creates a practical difficulty in that presently available 
pressure measuring techniques which are suitable for aircraft application 
are poor at measuring pressure differentials of less than 1/1200 
atmospheres which corresponds to an airspeed of 20 knots. 

Moreover, rotor downwash inevitably enters the pitot probe, which causes a fast 
error or enters the static port resulting in a slow error. Flow patterns developing around 
the aircraft during sideward and rearward flight may bring erroneous airspeed indications 
[Ref 13], 

In low speed flight, required engine power increases due to the difficulty of 
maneuvers. Due to these high engine and tail rotor anti-torque requirements, extra 
attention must be paid to directional control margins. Moreover, vibrator loads can occur 
in some maneuvers, which can result in fatigue damage accumulation in flight critical 
components [Ref. 4]. The pilot is often required to fly this type of maneuver, and ground 
references, if available, are mostly used instead of instruments. However, in a combat 
environment, ground references might not be available which increases the need for 
accurate measurement systems. 

Many developments have been completed since the 1950’s to increase the 
accuracy of measurement systems. One study suggested moving the probes above the 
rotor hub in order to protect the pitot system from the down-wash effect. Another study 
suggested using a swivel device mounted above the rotor hub that can measure true 
airspeed magnitude and direction [Ref 3]. In this design, two venturi tubes were mounted 
on opposite ends of a rotating arm to provide a differential measurement between the two 
sensors. These sensors were used to calculate airspeed and wind direction [Ref 4]. 
However, such a system requires a slip ring assembly or a similar means of transferring 
the data from the rotating system to the body of the helicopter. One other approach is 
based on the study of wake under the rotor and using a sensor mounted under the rotor to 
determine the airspeed of the helicopter. In this approach, a 360*^ rotating pitot-static 
probe was used to measure the true airspeed and wind direction [Ref 4]. 

In the late 1970s and early 1980s, Faulkner and Buchner proposed a study based 
on the idea that the measurement of the control position can be used to estimate airspeed 
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when airspeed substantially affects the control trim position. These researches built a 
model using simplified thrust and flapping equations to obtain longitudinal and lateral 
velocity components. However, their analysis requires knowledge of the main rotor flap 
angle, which is difficult to measure [Ref 8]. 

Although some progress has been made because of these studies, none of them 
have received worldwide acceptance due to system complexity, reliability, and economic 
and aerodynamic issues. 

McCool and Haas describe some other efforts to determine the airspeed 
analytically [Ref 3]. One of these efforts is a neural network-based approach, which is 
the core of this thesis. In that study, first the flight parameters, such as rotor RPM, cyclic 
position, etc. were obtained during test flights of HH-60J helicopter and CH-46E 
helicopter and then these parameters were used as inputs into a NN model to predict 
airspeed and wind direction. By using fuselage parameters, the problem of transmitting a 
signal from the rotating system to the fuselage is eliminated [Ref. 3]. The results obtained 
from this study proved that NN could be used to predict airspeed with reasonable 
accuracy. 

A NN implementation is chosen for several reasons. First, it is easy to use. Input 
parameters can be chosen as quantities, which are commonly measured on a flight data 
recorder. Second, the combined influence of inputs can be investigated by using multiple 
hidden layers. Third, a more flexible and empirical estimation can be obtained, as new 
data becomes available to retrain the network [Ref 16]. Based on the aforementioned 
benefits, McCool and Haas showed that the neural network approach provides a 
mechanically simple and inexpensive alternative to current low airspeed measurement 
technology [Ref 3]. 

The primary objective of this thesis is to determine the airspeed of an OH-6 A 
helicopter using an NN implementation with the input data provided from the helicopter 
simulator model FLIGHTLAB. Our intention is to increase the accuracy of the 
measurement systems in a cost-effective way by using the NN and a simulator model. 
Since simulator data is used as inputs to the system, first UH-60A data were obtained 
from the simulator and used as inputs to the NN. HH-60J and UH-60A helicopters are 
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very similar and they are both variant of Skorsky basie Hawk series helieopters. 
Therefore a eomparison of HH-60J helieopter and UH-60A helieopter results would 
provide the verifieation of using simulator data for the NN model. Onee the model is 
validated for the UH-60A helieopter, the next step is to prediet the speed of the OH-6A 
helieopter using the same approaeh. 
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II. NEURAL NETWORKS 


A. INTRODUCTION TO NEURAL NETWORKS 

The human brain is the world’s most complex computing device. Thus, the 
brain’s powerful thinking, remembering, learning and problem solving capabilities have 
been studied and modeled for digital computing. The brain is so complex that scientists 
have only just begun to understand how it works. Although little is known about the brain 
and the neural system, modeling the functionality of the brain in a very fundamental 
manner leads to the creation of neural networks. 

The neuron is the fundamental cellular unit of the nervous system and, in 
particular, the brain [Ref. 6]. The neuron functions like a microprocessing unit, which 
receives and combines signals from many other neurons. The brain consists of 10^^ 
number of neurons. They are connected to each other by approximately lO"^ connections 
per element or 60 trillion connections total [Ref. 1]. The input path of a cell body is 
called “dendrites”. The output path is called “axon”. The axon of a neuron splits up and 
connects to the dendrites of other neurons through a connection called a “synapse”. The 
cell body, which is also called “soma”, sums incoming signals when sufficient inputs are 
received. It is generally thought that all functions are stored in the neurons and in the 
connections between them. Learning occurs when new connections are established or 
existing connections are modified. The structure of a neuron is depicted in Figure 1 [Ref. 
10 ]. 

In an artificial neural network (NN), the same principles are used to simulate and 
capture some of the power of the brain and the neural system. In a NN, the unit 
corresponding to a neuron is called the “processing element (PE)”. A PE has many input 
paths and combines the values of these input paths. The combined input is then modified 
by a transfer function. 
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Figure 1. Structure of a Neuron [From Ref. 10]. 

The transfer function is also referred to as an “activation function”. This transfer 
function can be either a threshold function or a continuous function of the combined input 
[Ref 6]. When the transfer function is a threshold function, as when the combined 
activity level of PE reaches a certain level, the information passes or otherwise is filtered. 
The output path of the transfer function is passed directly to the output path of PE. The 
PE output is connected to the other PE input paths through connection weights, as in a 
neural system. 

The NN is trained to perform a particular function by adjusting the values of 
strength of connections (weights) between elements. The procedure to modify weights 
and biases of a network is called a learning rule, where the bias corresponds to a weight 
with a constant input. Note that, biases are sometimes added to add more flexibility to the 
network configuration. There are two types of learning; 

1) Supervised Teaming, where sets of inputs and desired outputs (target values) 
are presented to the network, and the NN configures itself to achieve desired input/output 
mapping; 2) Unsupervised Teaming, where only inputs are shown to the network and NN 
organizes itself internally so that each PE responds strongly to different sets of inputs 

In this thesis, only supervised learning schemes are considered. In that type of 

learning, a NN generates its own mles by learning from shown examples. Teaming is 

achieved through a learning mle that makes necessary modifications to weights and 

8 



biases in response to network output and target values. Figure 2 shows how such a NN 
system works. 



A NN consists of many PEs grouped together called layers. A layer is defined as a 
group of neurons having connections to the same inputs and sending outputs to the same 
destinations [Ref 2]. Figure 3 depicts the layer structure of a NN. 



Figure 3. Fayer structure of NN [From Ref. 10]. 

A general NN usually consists of an input layer, one or two hidden layers, and 
one output layer. The typical structure of a layer with one neuron is presented in Figure 4. 
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Figure 4. Structure of a layer [from Ref.2]. 


A layer consists of the following elements; 

• An input vector (p), 

• Weights (w) represent connection strength, 

• Summation (Z), 

• Bias vector (b) represents a column vector of bias values for a layer of 
neurons, 

• Transfer function (/), 

• Output vector (a) is the output of the network. 

The architecture of a NN depends on the number of layers the network possesses, 
the number of neurons in each layer, each layer’s transfer function, and how the layers 
are connected to each other. Note that there is no unique architecture for any given 
problem, and the best architecture heavily depends on the data presented to the network, 
the number of neurons for each layer selected by the user, and the selected activation 
function. Most practical NNs have just two or three layers. It is certain that the greater 
number of neurons in the hidden layer, the more powerful the network may potentially 
be. Note that, however, adding more neurons makes the network more complex and 
complexity should be minimized since it increases training time. Furthermore, the NN 
may memorize the training pattern set, but not perform well on data outside of the 
training set when the ratio of neurons to training patterns is too large. This problem is 
referred to as overspecialization or overfitting. 
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B, NEURAL NETWORKS BACKGROUND 

Although background studies of NN, extend baek to the late 19**' and early 20* 
eenturies, modern NN implementations first appeared in 1943 with the works of Warren 
MeCulloeh and Walter Pitts. Their researehes eonsidered the coneept of artifieial 
neurons, whieh have the eapability to eompute arithmetie or logieal funetions. MeCulloeh 
and Pitts published watershed paper entitled “A Logieal Caleulus of Ideas Imminent in 
Nervous Aetivity” [Ref. 6]. 

In the late 1950s, Frank Rosenblatt proposed the pereeptron network and the 
assoeiated learning rule. In 1957, Rosenblatt published the first major researeh projeet in 
neural eomputing whieh ineluded the development of the pereeptron element. The 
pereeptron is a pattern elassifieation system, whieh could identify both abstraet and 
geometrie patterns. In addition, the pereeptron ean make limited generalizations and ean 
properly eategorize patterns despite noise in the input [Ref 5]. This study showed the 
first praetieal applieation of NN by demonstrating how NNs ean perform pattern 
reeognhion. In 1959, Bernard Widraw and Tedd Hoff proposed the Adaline (Adaptive 
Linear Element), based on simple neuron-like elements and used it to train adaptive linear 
networks. The Adaline and the two-layer Madaline version were used for a variety of 
applieations ineluding speeeh reeognhion, eharaeter reeognhion, weather predietion, and 
adaptive eontrol. Widraw used the adaptive linear element algorithm to develop adaptive 
filters that eliminate phone line eehoes, in the first real life NN applieation [Ref. 6]. 

In the mid-1960s, Marvin Minsky and Seymour Papert, eonsidered the NN 
potential limitations. They showed that these already known networks eould handle 
linearly separable problems only and were usually not appropriate for real life 
applieations. As a result, NN researeh faded for a while. 

In the 1970s, Kohonen, Grossberg and Anderson proposed the Kohonen Network 
and the Self-organizing Network. Kohonen introdueed the eoneept of the eompetitive 
learning rule in whieh PEs eompete to respond to an input stimulus and the winner adapts 
itself to respond more strongly to that stimulus. This type eonstitutes an unsupervised 
learning proeess and the internal organization of the network is governed only by input 
stimuli [Ref. 5]. Grossberg’s eontribution was a wealth of researeh towards the design 
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and construction of neural model as he used neurological data to build neural computing 
models. Anderson developed a linear model, called a linear assoeiator. Which is based on 
models of memory storage, retrieval and recognition. In addition, Anderson improved the 
model by combining it with a nonlinear post-processing algorithm, whieh is used to elean 
up spurious responses. This model is ealled Brain-State-in-a -Box [Ref 5]. 

In the I980’s, NN became popular again with the back-propagation algorithm for 
training multilayer pereeptron networks. The concept of back-propagation algorithm was 
presented by several researchers, such as David Parker, Yaun LeCun, David Rumelhart, 
James McClelland, and Geoffrey Hinton. While a pereeptron network is only eapable of 
solving linear problems, back-propagation network ean solve more eomplex nonlinear 
problems. This significant capability made the baek-propagation networks the most 
widely used networks. 

In 1982, John Hopfield presented a paper describing his neural computing system 
ealled the “crossbar associative network” or known as the “Hopfield Model”. This model 
represented a neuron operation as a thresholding operation and illustrated memory as 
information stored in the interconneetions between neuron units. He also illustrated and 
modeled the brain’s ability to call up responses from many locations in response to a 
stimulus. Thus, this model represents how a NN associates information from many 
storage sites for a given input [Ref. 6]. In the 1980s, the Bi-direetional Assoeiative 
Memory (BAM) network, Boltzman Maehine, the General Regression NN, and the 
Learning Veetor Quantization Network were developed, in addition to the back- 
propagation and Hopfield models. 

Although the coneept of NN has been around for about 50 to 60 years, most 
applications have appeared in the last fifteen years and the field is still developing very 
rapidly. NNs can be found in many fields ranging from aerospace to medicine, banking 
and robotics. Given the work done and range of applications, NN will most likely be a 
permanent fixture not only as a solution to everyday problems but also as a tool to be 
used in appropriate situations. It is eertain that the more the structure of the brain is 
understood, the more advances there will be in NN. 
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C. INTRODUCTION TO THE BACK-PROPAGATION NEURAL 
NETWORK 


Although in 1974 Paul Werbos described a similar algorithm to train multilayer 
network, Back-propagation (BP) was invented independently by David Parker, Yaun 
LeCun, David Rumelhart, Geoffrey Hinton and Ronald William in the mid 1980s. 

Back-propagation is a very popular NN learning algorithm to train multilayer 
networks with differentiable transfer functions to perform function approximation, 
pattern association and classification. The algorithm is named for how it handles errors in 
the network. The derivatives of the network error with respect to the network bias and 
weights are computed starting from the very last layer up to the first layer in order to 
modify the weights and biases of the network. The process starts from the output layer 
and goes back to the first layer and is therefore called the “Back-Propagation Network”. 

Typically, back-propagation network has an input layer, one or two hidden layers 
and one output layer. Although there is no limit on the number of hidden layers, generally 
one or two hidden layered networks are selected due to the complexity of the resulting 
system. An example of a typical network is shown in Figure 5. 
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Figure 5. A Typical Back-Propagation Network [From Ref 2]. 

12 3 

A BPNN is typically represented with the following notation; R-S -S -S , where 

R is the number of inputs and S' is the number of neurons at layer i. In Figure 5, the input 

vector to the network is shown by P\. Transfer function is represented by / and 1) 

represents bias. Depending on the selected transfer function, the BPNN can be used to 

13 











































solve linear or nonlinear problems. In order to solve a nonlinear problem, a nonlinear 
transfer function should be selected. 

The first step in BPNN is to transfer the input forward through the network. The 
output of one layer becomes the input to the following layer. Thus, the output of a 
network can be defined as 


=f^{w^ +b^) , ( 2 . 1 ) 

where M = 1,.,R and R is the number of layers in the network. In Figure 5, for the 

first layer, is a*^^^ which is the input of the network, and a*^^^ is the output of the last 
layer. 

The second step is to propagate the errors, or in other words, sensitivities 
backward through the network from the last layer all the way to the first layer. The 
BPNN algorithm is a generalization of the LMS algorithm, which uses the mean square 
error as a performance index. As each input is applied to the network, the output, a, is 
compared with the associated target value, t. The algorithm then adjusts the network 
parameters in order to minimize the mean square error [Ref. 1]. 

Thus, given the error is defined as; 

e=\t-d\ , (2.2) 

and the expectation of mean square error is defined as 

fix) = Eiie^e]) = E[it - af {t - a)], (2.3) 

the expectation of the squared error at iteration k becomes [Ref 1] 

fix) = Eiitik) - aik)y itik) - aik ))-\. (2.4) 


The algorithm adjusts weights and biases as follows [Ref 1] 

df 


Wij (k+l) = Wi. (k)-« 




h.^(k+l) = b“(k)-a 


db!^ 


(2.5) 


14 




where a is the learning rate. This shows that the weights at any given iteration are equal 
to the weights at the previous iteration adjusted by some fraetion, a, of the sensitivity of 
the error to that weight. In other words, the weights at eaeh iteration are adjusted in a way 
that reduees the error at the previous iteration. 

The seleetion of a is usually done by trial and error. Too large an a often leads to 
divergence of the learning algorithm while too small an a results in a slow learning 
process [Ref 9]. 

Note that the above equations include partial derivatives. Since the error is an 
indirect function of the weights in the hidden layer, a chain rule is used to compute these 
partial derivatives. These partial derivatives are called sensitivities and defined as: 


M 


dnf 


( 2 . 6 ) 


In the above equation Si is called the sensitivity of f to changes at layer M in the i* 
element of the net input. 


The next step is to propagate the sensitivities backward 


. M-l 

n^-\w^Ys^, forM = R,....,l . 


(2.7) 


The final step is to adjust weight and biases with respect to these sensitivities, which 
leads to 


{k + \) = w" {k) - as^ , 


( 2 . 8 ) 


b^(k + \) = b^(k)-as^ . 

Although the BP network is a very powerful technique, it usually requires a long 
training time. However, using some tricks and heuristics can often help to obtain an 
efficient BPNN implementation [Ref. 9]. 

In summary, the BPNN is a commonly used technique for solving nonlinear 
estimation problems, provided that a nonlinear transfer function is selected. The BP 
algorithm is based on a LMS algorithm that minimizes the squared error, where the chain 
rule is used to compute the derivatives of the error with respect to weights and biases in 
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the hidden layers. The methodology in the BP network is a three-step proeess. First, the 
input is propagated forward through the network and then the sensitivities are eomputed. 
Seeond, starting from the last layer, these sensitivities are propagated baekward through 
the network. Finally, weights and biases are updated using these sensitivities. 
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III. DEVELOPMENT OF NEURAL NETWORK MODEL 


A, NEURALWORKS PROFESSIONAL II/PLUS SOFTWARE 

The model built for this analysis is derived from the NeuralWorks simulation 
software program presented by NeuralWare, Ine. NeuralWorks is a user-friendly program 
that makes it possible to not only seleet network parameters easily and quiekly but also to 
present network results effeetively. Many NN algorithms sueh as BPNN, Radial Basis 
Funetion NN, and LVQ, are ineluded in the software. Onee the algorithm type is ehosen, 
the next step is to define the network arehiteeture, whieh ineludes speeifying the number 
of layers and the number of PEs assoeiated with eaeh layer, etc. Several types of learning 
rules and PE transfer functions are embedded in the program. NeuralWorks also allows 
the user to select learning rates and momentum terms. In addition, the extensive and 
powerful instrumentation and diagnostic package allows the user not only to monitor and 
adjust network parameters but also to display weights, errors, classification rates and 
confusion matrices in graphical formats. Moreover, users can display the networks 
through network or Hinton diagrams. These specifications make NeuralWorks 
Professional II/Plus useful in designing, building, training, testing and deploying neural 
networks to solve complex, real-world problems [Ref 5]. 

The following picture shows a typical window of the Neural Works Professional 
II/Plus software; 



Eigure 6. Typical Window of Neural Works Professional II/Plus Software. 
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B, FLIGHTLAB SIMULATOR 

Flight simulation tools have been used extensively not only for training and 
evaluation of airerew members but also for the design and analysis of aerospaee vehieles. 
The need to use flight simulation systems instead of real aerospaee vehieles depends on 
the following reasons: 

• The eomplexity of aireraft systems, 

• High operating eosts of airerafts, 

• The limitations of operating environment of airerafts, 

• Teehnologieal improvements in flight simulators. 

In addition, flight simulators provide safe and effeetive eonditions for training 
purposes as instruetion, demonstration and praetiee of eertain maneuvers and proeedures 
that eannot be done during real flight eonditions may easily be performed in a simulated 
environment. Finally, longer training periods ean be tolerated due to the low operating 
eost of simulators. Therefore, modem training proeedures beneht from simulation tools 
extensively. 

Simulators also play an important role in engineering proeesses sueh as design 
and analysis of aireraft. It is eertain that the dynamies of rotary-wing aireraft is mueh 
more eomplex than those of fixed-wing aireraft. Consequently, helieopter simulation 
systems require high eomputational power to inelude a eomplex set of eomputer 
programs, very powerful and expensive eomputer systems and full motion based 
simulation deviees. Reeent improvements in eomputer and simulation teehnology make 
simulators able to produee data as eorreet as real helieopters. These improvements make 
it possible to make use of simulators for design, test and analysis purposes. 

FLIGHTLAB is a eommereial-off-the-shelf type of software produet developed 
by Advaneed Rotoreraft Teehnologies, Ine (ART) for the development and operation of 
flight vehiele dynamies models in simulation applioations espeeiahy for helieopters. 
FLIGHTLAB is a tool that allows a user to build eaeh seetion of a model separately and 
then combine the pieces under the same framework similar to building a model by using 
a finite element approach [Ref 14]. This approach is very useful for engineering studies. 
In order to provide a desktop pilot interface for FLIGHTLAB models, ART has 
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developed PilotStation, a turnkey exeeutive for real time simulation operations that 
eouples image generation and pilot eontrol inputs with the FLIGHTLAB flight dynamies 
model to provide a low eost, pilot-in-the-loop simulation. PilotStation ean also be used 
with FLIGHTLAB eode generated models to provide a low eost, real time simulation 
eapability using PCs [Ref 7]. PilotStation, whieh ean be run on LTNIX or LINUX 
operated eomputers, proeesses the high fidelity rotoreraft model and generates the eoekpit 
gauges and the window displays. 

FLIGHTLAB uses a simulation language ealled Seope. This language is an 
interpretive language that uses the industry standard MATLAB syntax, eoupled with new 
language eonstruets and a eombination of C and Fortran eomputer languages to make the 
building, testing and solving of nonlinear dynamie rotoreraft simulation models easy 
[Ref 7]. 

Another feature of the FLIGHTLAB Simulation System is Xanalysis. This tool 
makes it possible to design, test and analyze rotoreraft models by allowing the user to 
modify model parameters and perform a wide range of analyses on design alternatives. A 
set of predefined test seenarios earry speeifie rotoreraft analyses sueh as performanee, 
stability and eontrol, handling qualities and aerodynamie and struetural loads. In addition 
to these features, the simulation ean be automatieally eonfigured to interaet with the test 
vehiele’s eonfiguration and test eonditions while using the time history of the test flight’s 
eontrol inputs to implement the simulation [Ref. 7]. These features are the required 
eharaeteristies neeessary to implement the model built for this thesis and make use of 
data provided by a simulator instead of a real test flight. 

The following figure depiets the X-Analysis "Flight Test Utility" [Ref. 7]; 
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Figure 7. Xanalysis Window [From Ref. 7]. 

C. SELECTING DATA AND BUILDING THE MODEL 

In the early 1990s, Kelly M. MeCool and David J. Haas at Naval Surfaee Warfare 
Center, Bethesda, MD, developed a study of UH-60 helieopter airspeed estimations using 
NN in a low speed environment. McCool and Haas built and implemented a BPNN with 
two hidden layers and 25 processing elements in each layer. The data used as input for 
the NN were obtained from actual test flights performed at the Naval Air Warfare Center, 
Aircraft Division, Patuxent River, MD. [Ref 4]. Their study was used as a starting point 
for this thesis since similar analysis was performed but with different data sets. They used 
real flight test measurements as a source of data for their model. FLIGHTLAB simulator 
outputs were used for NN input for this thesis. In order to evaluate how well the 
simulator performs and how well the relationship between the simulator parameters and 
low airspeed is developed, the simulator was first run for the UH-60A helicopter to 
compare the results of the present study and the results of the actual test flights. The 
simulation data was analyzed using a similar network architecture and network 
parameters, to make a fair comparison. The goal of this approach is to determine whether 
using a simulator instead of a real test flight is feasible or not. 

The simulator was run for UH-60 A helicopter from hover to 50 knots with 5 knots 
intervals at various gross weights ranging from 16000 to 24000 lbs with 1000 lbs 
intervals and at various sideslip angles from 0 to 360*^ with 30° intervals to obtain the 
input data for the NN model. However, sideslip angles were varied from 300° to 60° for 
velocities 35 knots and above because of the limitations of the UH-60 A helicopter since 
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the helicopter cannot fly rearward or sideward when the speed is above 30 knots. The 
altimeter was set to 85 feet and AGL to obtain the results when the helicopter is out of the 
ground effect and at level flight. In a ground effect, analysis was performed by setting the 
altimeter to 20 feet AGL for the UH-60A helicopter. The wind was assumed to be zero 
for this analysis. Fifteen parameters were chosen as inputs for the neural network model 
based on the variables described above, and are shown in Table 1. 


MODEL VARIABLES FOR UH-60A HELICOPTER 

Airspeed 

0-50 Knots (5 Knots intervals) 

Gross weight 

16000-24000 lbs (1000 lbs intervals) 

Sideslip angle 

0- 360 degrees and 300 - 60 degrees (30 degrees intervals) 

Altimeter (AGL) 

20 ft (for in ground effect), 85 and (for out of ground effect) 

Wind Speed 

0 

Pressure Altitude 

Sea level 

NN MODEL INPUT PARAMETERS 

1. 

Cyclic position (Lateral) 

2. 

Cyclic position (Longitudinal) 

3. 

Collective position 

4. 

Pedal position 

5. 

Roll rate 

6. 

Pitch rate 

7. 

Y aw rate 

8. 

Pitch attitude 

9. 

Roll attitude 

10. 

Altimeter 

11. 

Climb rate 

12. 

Main Rotor Blade RPM 

13. 

Engine torque 

14. 

Gross weight 

15. 

Sideslip angle 


Table 1. Neural Network Input Parameters. 
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These fifteen parameters, whieh define the eondition of the related part of the 
helieopter at eertain settings, are the outputs of the FLIGHTLAB simulator. The data was 
split into two sets, one for training the network (training data set) and the other one for 
evaluating the network performanee (test data set). The first set used for training ineludes 
the data related to airspeeds at 0, 10, 20, 30, 40, 50 knots at all gross weights and sideslip 
angles. The seeond set used for testing was formed with the remaining input data. 

The type of network was seleeted, after setting up training and test data sets. A 
Baek-propagation approaeh was seleeted as it was used before with sueeess [Ref 3]. A 
Radial Basis Funetion Network (RBFN) was seleeted as an alternative network beeause it 
ean be used in most applieations where baek-propagation approaehes may be used. 
However the RBFN trains faster and leads to better deeision boundaries than a BPNN in 
many elassifieation and deeision problems. The learning rule used in RBFN is an 
unsupervised learning rule. The results of these two models are presented in Chapter IV. 

The next step in the BPNN seheme involved the seleetion of number of layers and 
neurons in eaeh layer. First the model in Referenee 3 was used as a starting point, sinee 
the first goal of this study was to eompare the results of the model implemented by real 
flight data and the model by simulator data. Thus, we seleeted two hidden layers and 25 
neurons in eaeh hidden layer, resulting in a 15-25-25-1 BPNN strueture. Several models, 
were also studied; 14-25-25-1, 15-18-25-1, 14-14-14-1, 14-14-12-1, and 14-14-10-1. 

Next, the learning eoeffieient and momentum term were determined. Learning 
eoeffieients eontrol the ehanges in size of weights and biases during learning. Setting an 
appropriate learning rate is signifieant beeause a small learning eoeffieient leads to very 
little learning, whieh inereases the training time. However, a large learning eoeffieient 
may eause the performanee index to diverge, meaning that no learning oeeurs. Therefore, 
a small learning rate is generally used to avoid divergenee. Next, the momentum term is 
seleeted. The momentum term helps to obtain faster learning when using a low learning 
rate. In this study, the appropriate values for these two terms were determined by trial and 
error. 

Finally, the learning rule and transfer funetion type were seleeted. There are 
several learning rules available in NeuralWorks, sueh as, the Extended Delta-Bar-Delta 
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(Ext. DBD), Normalized Cumulative Delta (NCD), Quickprop, Delta Bar Delta (DBD) 
and Delta Rule (DR). In this study, several types of learning rules were analyzed for the 
UH-60A helicopter and the results are presented in Chapter IV. Based on the results of 
the UH-60A helicopter model, only Ext. DBD and NCD learning rules were used for the 
analysis of the OH-6A helicopter. NCD is a learning rule, which is immune to changes in 
the epoch length, where an epoch is the number of sets of training data presented to the 
network between weight updates. The Ext. DBD rule takes momentum term into 
consideration. The transfer function is used to map the output of a neuron or a layer to its 
actual output [Ref. 2] and may be linear or nonlinear. The most commonly used transfer 
functions are depicted in Table 2 [Ref 1]. 


NAME 

INPUT/OUTPUT RELATION 

Hard Eimit 

a = 0 n < 0 


a = 1 n > 0 

Einear 

a = n 

Eog-Sigmoid 

1 

a =- 

1 + A" 

Hyperbolic Tangent Sigmoid 

n -n 

e -e 

a = - 

e +e 

Competitive 

a = \ neuron with max n 


a = 0 all other neurons 


Table 2. Commonly Used Transfer Eunctions. 


The sigmoid transfer function is commonly used in BPNN implementations 
because it is differentiable. NeuralWorks provides the hyperbolic tangent and the sine 
function as alternative functions. The hyperbolic tangent function is a bipolar version of 
the sigmoid function. The sigmoid function maps the output between 0 and 1 smoothly 
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whereas the hyperbolie tangent funetion plaees it between -1 and 1. In this work, the 
hyperbolie tangent transfer funetion (Tanh) is used. 

Onee the feasibility of this approaeh had been verified, the same type of 
arehiteetures and NN model parameters, were used for the OH-6A helieopter. These 
parameters were the same as those given in Table 1, exeept for the engine torque 
parameter. The only difference utilized for the second part of the analysis was the 
altimeter setting. The altimeter was set to 12 feet for in ground effect data and 100 feet 
for out of ground effect data for the OH-6A helicopter. Besides the altimeter setting, 
other parameters, such as gross weight and pressure and altitude, were adjusted for the 
OH-6A helicopter. Gross weight for OH-6A was defined ranging from 1500 lb. to 2500 
lb. with 100 lb. increments. Pressure altitude was set to 90 ft for sea level and 6000 ft for 
high-level altitude. For both pressure altitudes, in ground effect analyses were performed 
at 12 ft, while out of ground effect analyses were performed at 100 ft. 
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IV. ANALYSIS AND RESULTS 


The NN-based approach to determine helicopter low airspeed developed by Haas 
and McCool has proven to be quite promising. In their study, a BPNN configuration with 
the Extended Delta Bar Delta learning rule, 2 hidden layers and 25 PEs per layer was 
selected [Ref. 3]. The data used as input to the NN were obtained from actual test flights 
performed at the Naval Air Warfare Center, Aircraft Division, Patuxent River, MD. 
Based on the model they implemented, and using the methodology described in the 
previous section, in this thesis the NN model was first developed for the UH-60A 
helicopter by using simulator data. The first model architecture developed for this study 
was approximately the same model as that described in Ref 3, since most of the network 
parameters such as number of layers and number of neurons are determined by trial and 
error. The first network selected was a 14-25-25-1 BPNN with the Ext. DBD learning 
rule, and did not include the sideslip angle parameter as input to the network. Various 
architectures were subsequently tried with different learning rules. The momentum 
function was varied from 0.2 to 0.6 with the best results obtained at 0.4. Therefore, all 
OH-6 A analyses used a momentum function to 0.4. The NN output error airspeed was 
measured in knots at 1 a and in terms of the root mean square (RMS) error, as reported 
by the NeuralWare software. The number of learning iterations was set to 50000 since the 
network stabilized after approximately 20000 iterations for most of the implementations. 

A helicopter flying close to the ground requires less power than when it is flying 
far from the ground. The proximity of the ground to the rotor disk constrains the rotor 
wake and reduces the induced velocity at the rotor, which causes a reduction in the power 
required for a given thrust. This is called ground effect. A helicopter can hover in ground 
effect at a higher gross weight or altitude than when it is out of ground effect. However, 
in forward flight, the effect of the ground diminishes with the forward speed. Data sets of 
both helicopters were obtained by running the simulator for two flight conditions: in- 
ground effect (IGE) and out-of-ground effect (OGE). These data sets are called single 
condition data sets. After training and testing with the single condition data sets 
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separately, the sets were eombined for a baseline data set. The networks were retrained 
and tested with the eombined data set. 

The OGE single eondition data set has a dimension of [1144x14] for the UH-60A. 
The IGE single eondition data set has a dimension of [1144x15]. The differenee is due to 
the faet that the EiH-60A OGE data did not inelude the sideslip angle. Sideslip angle was 
ineluded for the EIH-60A IGE and eombined data sets. Eor the OH-6A helieopter type, 
the sideslip angle was ineluded in all data sets, but the engine torque was not available, so 
all OH-6 A single eondition data sets have dimensions of [1144x14]. The single eondition 
data sets were split into training and testing sets. Eaeh training set is a subset of its single 
eondition data set with a dimension of [638x14] or [638x15]. Eaeh testing set is made up 
of the remaining elements of the single eondition set and has a dimension of [506x14] or 
[506x15]. Various networks were trained separately using these single eondition data 
sets. Eater, the single eondition data sets were eombined forming a baseline data set with 
a dimension of [2288x15] for the EiH-60A and [2288x14] for the OH-6A. A baseline 
training set was formed as a subset of the eombined set with a dimension of [1276x15 or 
14]. The remaining elements of the baseline set formed the baseline testing set with 
dimensions of [1012x15 or 14]. The networks were retrained using these eombined 
training sets. Combined and single eondition, OGE and IGE, results were obtained again 
from eaeh network that was trained with these baseline data sets. Einally, NN model 
results were exported to MATLAB to develop the following tables and figures. 

Eaeh figure eonsists of four windows. The first window shows the relationship 
between gross weight and the predieted speed related to that partieular gross weight. The 
seeond window shows the range of the predieted speed for eaeh aetual speed. The third 
and fourth windows display the relationship between predieted speed and the sideslip 
angle. In the third window, the speed range is from 0 knots to 30 knots, whereas in the 
last window, the speed varies from 35 knots to 50 knots. In order to make the figures 
more readable and easier to distinguish the NN predictions from one another, adjacent 
speeds are illustrated with a different marker. NN outputs for 5 knots, 15 knots and 45 
knots are presented with circles while outputs for 15 knots and 35 knots are demonstrated 
with triangles. The structure of the NN, learning rule, helicopter type, RMS error and 

altitude information about that figure is provided in the figure label. 
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Each table corresponds to a different NN type and consists of five column s . The 
first column presents actual airspeed and target values for the NN model. The second 
column gives the mean value of the NN outputs related to each target speed. The third 
column shows airspeed errors at 1 a in terms of knots. This error shows the standard 
deviation (SD) around the mean and it is computed by the MATLAB function “std”. The 
next column demonstrates the error percentage at 1 a relative to the target speed. It is 
computed by the following formula: 


Percent Error = — * 100 , 
Vi 


(4.1) 


where Gi is the error SD related to that speed and Vi is the actual speed. The last column 
shows the absolute maximum error of the NN predicted speeds. 

Results obtained for analyzed architectures follow in the next section. 

A, ANALYSIS OF THE UH-60A MODEL 

1, Out of Ground Effect (OGE) Analysis 

This section shows the results of OGE analyses by setting the altimeter of the 
simulator to 85 feet AGE. The network was trained using this single condition data set 
obtained for this altitude. The results of OGE analysis using the combined data set are 
shown in the UH-60A Baseline Analysis section. 

a. 2-Hidden Layer BPNN 

(1) 14-25-25-1 Ext. DBD, Results are shown in Eigure 8 and in 

Table 3. The RMS error is 0.0593. Although the network performance on the training 
data is quite good, the test results are not close enough to the target speeds. Note that, the 
NN produced a maximum speed error of 3.51 and 3.86 knots at 5 knots and 45 knots 
respectively. A 0.7017 knots airspeed error SD was achieved by this setup. The 
percentage of error is worst at 5 knots, where it is 8.5%, while it is about 3 % for other 
speeds. 
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Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7017 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.3322 

0.4249 

8.4979% 

3.5125 

15 

16.1533 

0.4601 

3.0675% 

1.8307 

25 

23.8216 

0.5045 

2.0179% 

2.2057 

35 

36.2752 

1.0928 

3.1222% 

3.1262 

45 

46.2789 

1.4228 

3.1617% 

3.8625 


Table 3. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

Ext. DBD learning rule. 
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Figure 8. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

Ext. DBD learning rule. 


(2) 14-25-25-1 NCD. Results for this configuration are shown 

in Table 4 and in Figure 9. The RMS error is 0.0658 with the NCD learning rule. The 
maximum value of the predicted speed error is ±3.8584 knots. The speed prediction 

accuracy within ± 1 a is 0.7506 knots. The airspeed error SD at speeds 25, 35 and 45 
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knots is about 1 knot or 3 %. For speeds equal to 5 and 15 knots, the error percentage 
increases to 6 %. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7506 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.7932 

0.3221 

6.4425% 

2.7872 

15 

17.0481 

0.8209 

5.4726% 

3.0500 

25 

23.3043 

0.6470 

2.5878% 

2.5837 

35 

36.8160 

1.0996 

3.1417% 

3.8584 

45 

46.2666 

1.1273 

2.5051% 

3.0720 


Table 4. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

NCD learning rule. 
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Figure 9. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

NCD learning rule. 
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(3) 14-25-25-1 DBD. Figure 5 and Table 10 summarize the 

findings for this configuration. Results show the RMS error for this learning rule is 
0.0721. This model predicts the speed very accurately when the actual speed is 45 knots 
with a mean value of 45.19 knots, corresponding to a speed error prediction rate of 


1.42%. However, when the speed is 15 knots and 35 knots the maximum error in the 
prediction increases significantly to 4.9681 knots, while the error SD goes up to 1.1142 
knots. The total airspeed error for the DBD configuration is equal to 0.7544 kn ots. 
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Figure 10. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 


DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7544 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.8389 

0.1628 

3.2559% 

0.4681 

15 

11.3654 

0.5467 

3.4647% 

4.5417 

25 

24.5698 

1.0986 

4.3944% 

2.2028 

35 

38.0388 

1.1142 

3.1834% 

4.9681 

45 

45.1962 

0.6442 

1.4227% 

1.2996 


Table 5. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 


DBD learning rule. 
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(4) 14-25-25-1 QUICKPROP. Figure 11 and Table 6 present 

results obtained for this configuration. The RMS error for the quickprop learning rule is 
0.0711. The maximum predicted speed absolute value error is 4.6189 knots. Note that the 
results are very close to those obtained for the DBD configuration. Although the mean 


airspeed value at 15 knots is significantly lower than that of the target speed, the error SD 


is the largest at 25 knots. The overall network airspeed error SD is 0.7326 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7326 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.5007 

0.1677 

3.3546% 

0.8191 

15 

11.4068 

0.6067 

4.0444% 

4.6189 

25 

24.8524 

1.1168 

4.4673% 

2.3644 

35 

37.9663 

0.8684 

2.4811% 

4.4915 

45 

44.6973 

0.5294 

1.1764 

1.1524 


Table 6. 


Results 


br the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 


Quickprop learning rule. 


Predicted Speed vs Grossweight 
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Figure 11. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

Quickprop learning rule. 
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b. 1-Hidden Layer BPNN 

Several one hidden layer BPNNs were implemented; the configurations 
considered were 14-8-1, 14-10-1, 14-12-1, 14-20-1, 14-22-1, 14-25-1, and 14-50-1. 
Results for best performing networks are shown below. 

(1) 14-20-1 NCD. Table 7 and Figure 12 summarize the 

findings for this configuration. The 14-20-1 NCD model performs quite well for the 
velocity of 25 knots with the error SD of 0.5 knots and percent error of 2 %. For all other 
velocities, the error percentage is 4 %. The RMS error is found to be 0.0407. The 
maximum error is 3.3 knots at the speed of 5 knots and 35 knots. However, the accuracy 
of prediction is better than many 2-hidden layer BPNN. The overall error SD is ± 0.8394 
knots. 
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Figure 12. Results for the UH-60A helicopter at 85 ft.; network configuration 14-20-1; NCD 

learning rule. 
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Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.8394 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.6992 

0.2008 

4.0170% 

3.3712 

15 

13.7729 

0.6762 

4.5079% 

1.6789 

25 

25.0478 

0.5373 

2.1494% 

1.2362 

35 

35.6197 

1.5618 

4.4623% 

3.3712 

45 

45.4842 

1.5915 

3.5366% 

2.8769 


Table 7. Results for the UH-60A helieopter at 85 ft.; network eonfiguration 14-20-1; NCD 

learning rule. 

(2) 14-20-1 Ext. DBD, Table 8 and Figure 13 summarize the 

findings obtained for this eonfiguration. Results show the RMS error is 0.0657 and that 
the error SD increases along with the speed. The maximum error is 3.7970 knots for a 35 
knots speed. The error SD is also the highest at this speed. However, results show a 
significant degradation of the airspeed prediction in % at 5 knots. The overall error SD is 
0.7002, which is quite better than that obtained with the other networks. 
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Figure 13. Results for the UH-60A helicopter at 85 ft.; network configuration 14-20-1; Ext. 


DBD learning rule. 
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Aetual 

Airspeed (kts) 

NEURAF NI 
Total 

iTWORK RESUFTS 

1 SD = 0.7002 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Pereent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.8849 

0.4173 

8.3461% 

2.9556 

15 

17.4373 

0.7170 

4.7798% 

3.2811 

25 

23.7263 

0.6597 

2.6389% 

2.3582 

35 

36.9723 

1.0187 

2.9105% 

3.7970 

45 

46.1386 

0.9508 

2.1129% 

2.6505 


Table 8. Results for the UH-60A helieopter at 85 ft.; network eonfiguration 14-20-1; Ext. 

DBD learning rule. 

2. RBFN Networks 

Figure 14 and Table 9 present findings obtained for the 14-200-1 RBFN 
eonfiguration. Several RBFN network eonfigurations were implemented. The best result 
was obtained with the 14-200-1 NCD eonfiguration; where the RMS error is 0.0733. 
However, results also show a larger veloeity varianee and maximum error than those 
obtained with BPNN implementations. In addition, the maximum speed error is 9.2719 
knots, whieh is mueh larger than that observed with BPNN implementations. The error 
SD is around 1 knot at all speeds, exeept at 45 knots where it signifieantly inereases to 
about 4 knots. The overall error SD is 1.634. Therefore, results showed that BPNN 
implementations are better suited than RBFN eonfigurations for this problem. 


Aetual 

Airspeed (kts) 

NEURAF NETWORK RESUFTS 

Total SD= 1.6340 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Pereent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.1368 

1.0915 

21.8301% 

5.4400 

15 

14.3318 

1.0551 

7.0338% 

2.7856 

25 

25.6402 

1.1621 

4.6485% 

2.9813 

35 

35.9359 

1.4807 

4.2307% 

3.8438 

45 

44.0553 

3.8565 

8.5699% 

9.2719 


Table 9. Results for the UH-60A helieopter at 85 ft.; network eonfiguration 14-200-1; 

RBFN NCD learning rule. 
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Figure 14. Results for the UH-60A helieopter at 85 ft.; network configuration 14-200-1; 

RBFN NCD learning rule. 


3. Pruned BPNN 

Figure 15 and Table 10 present the findings for this configuration. The BPNN was 
implemented with a pruning capability embedded in the NeuralWare software to increase 
the prediction efficiency, where pruning is used to remove connections from the network 
when their contributions are very small [Ref 6]. 

In light of the previous analyses, the 14-25-25-1 BPNN configuration with the 
NCD learning rule was chosen. In this scheme, data were presented to the network as 
before but the network performance was monitored against the previous ones during the 
training process after each 1000 iterations. The RMS error was compared at each 
“checkpoint”, and PE contributions checked. At that point the training stopped when the 
RMS was larger than at the last check, otherwise, it continued for the next 1000 
iterations. During this process, PEs with contributions smaller than a given tolerance 
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were disabled, which reduces network complexity and increases generalization efficiency 
[Ref. 4], 


In this analysis, the process went up to 46000 iterations and 7 PEs in the first 
hidden layer were disabled. Results show that the RMS error is 0.0355 which is best 
among all architectures investigated. The maximum predicted speed error is ±3.35 knots, 
which occurs when the actual speed is 35 kn ots. Errors for low velocities are below ±1.8 
knots. The error percentage for 5, 15 and 35 knots is about 4 % and about 2.5 % for other 
speeds. Results show that the predicted speed mean values are very close to the target 
values when compared with previous networks, finally, the overall error SD is 0.7167, 
which is quite good. 
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figure 15. Results for the EIH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

NCD (Pruned) learning rule. 


36 



















Actual 

Airspeed (kts) 

NEURAE NETWORK RESUETS 

Total SD = 0.7167 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.6789 

0.2113 

4.2253% 

1.7089 

15 

14.5970 

0.6322 

4.2147% 

1.6131 

25 

24.9281 

0.6610 

2.6440% 

1.2897 

35 

36.0136 

1.3438 

3.8395% 

3.3496 

45 

46.0334 

1.0039 

2.2309% 

2.6618 


Table 10. Results for the UH-60A helicopter at 85 ft.; network configuration 14-25-25-1; 

NCD (Pruned) learning rule 

The above results showed that the networks with NCD and Ext. DBD learning 
rules and pruned network have better performance than those obtained with other 
networks. Therefore, all other analyses, UH-60A in-ground effect, baseline data and OH- 
6A analyses, were performed using these configurations only. 

4, In-Ground Effect Analysis 

In-ground effect analysis of the UH-60A helicopter is performed when the 
pressure altitude is at sea level and the altimeter above ground is at 20 feet. The sideslip 
angle parameter was added to the NN inputs. Three network architectures were used for 
the in-ground effect analysis; 15-25-25-1 BPNN with NCD and Ext. DBD learning rules, 
and 15-18-25-1 BPNN with the NCD learning rule. 

a. 15-25-25-1 BPNN NCD 

Eigure 16 and Table 11 summarize the findings for this configuration. The 
RMS error is 0.0374 and the maximum error is 3 knots (at 35 knots). The error SD is also 
the highest at this speed. The network performance at low speeds and at 45 knots is quite 
good, with the error SD of less than 1.2, but the error SD goes up to 1.8 knots for a 35 
knots speed, which corresponds to about 5 %. The overall airspeed error SD is 0.8469 
knots. 
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Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.8469 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.1567 

0.2334 

4.6686% 

1.3234 

15 

14.0652 

0.5405 

3.6033% 

2.1345 

25 

25.8750 

0.6741 

2.6963% 

2.5196 

35 

35.5357 

1.8206 

5.2017% 

3.0036 

45 

45.2548 

1.2123 

2.6941% 

2.2981 


Table 11. 


Figure 16. 


Results for the UH-60A helicopter at 20 ft.; network configuration 15-25-25-1; 

NCD learning rule. 


Predicted Speed vs Grossweight 


Predicted Speed vs Actual Speed 


45 

^35 

•o 

V 

<v 

- 25 

•o 

4* 

tJ 

I 15 

a. 

5 

0 


30 
? 25 


mmm 

^45 

2 

35 

1 


& 

i25 

$ 


I 15 

a. 


000000000 

5 

0 



15 

5 

0 


1.6 1.8 2 2.2 2.4 

Gross Weight (lb*10000) 

Slow Speed vs Angle 


qOOqOOOo^ OOO 
< < <] <1 < <10 <] <1 < < 0 
OOOOOOOOOOOO 


0 5 


45 


15 25 35 

AHual Speed (HI) 

Fast Speed vs Angle 


60 120 180 240 300 360 

Side Slip Aigle (Deg) 

15-25-25-1 BPNN NCD UH-60A RMS = 



Side Slip Aigle (Deg) 

0.0374 ALT = 20 ft Be a Level! 


Results for the UH-60A helicopter at 20 ft.; network configuration 15-25-25-1; 

NCD learning rule. 

b. 15-25-25-1 BPNN Ext DBD 


Figure 17 and in Table 12 present the results for this configuration. The 
RMS error obtained is 0.0637. Note that the maximum error occurs when the helicopter is 
moving at 45 knots with a -60° angle, unlike other networks. Moreover, the maximum 
error, 4.81 knots, is higher than that obtained in most networks. Apparently, the error SD 
is also larger at 45 knots. Although the error SD at low speeds is less than 0.88 knots, the 
error percentage is about 6%. Overall, the error SD is 0.8308 knots. 
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Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.8308 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.5383 

0.3051 

6.1027% 

3.1938 

15 

16.7466 

0.8875 

5.9166% 

3.7167 

25 

23.4740 

0.5188 

2.0750% 

2.1472 

35 

35.6272 

0.9530 

2.7229% 

2.1010 

45 

46.5186 

1.6658 

3.7018% 

4.8112 


Table 12. Results for the UH-60A helieopter at 20 ft.; network configuration 15-25-25-1; 

Ext. DBD learning rule. 
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Figure 17. Results for the UH-60A helicopter at 20 ft.; network configuration 15-25-25-1; 

Ext. DBD learning rule, 
c. 15-25-25-1 BPNNNCD (Pruned) 

Table 13 and Figure 18 present the results for this configuration. Seven 
PEs were disabled in the first hidden layer. The RMS error is about 0.0724. The 
maximum error is 5.2 knots at 35 knots. The error SD is also larger at this speed. The 
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overall error SD is 0.8884 knots which is slightly larger than the NCD 15-25-25-1 
network. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.8884 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

4.7761 

0.2380 

4.7603% 

0.7986 

15 

11.4320 

0.4591 

3.0609% 

4.8662 

25 

24.8531 

1.3065 

5.2260% 

2.8201 

35 

37.9958 

1.4152 

4.0435% 

5.2195 

45 

45.4203 

0.7768 

1.7263% 

1.9679 


Table 13. Results for the UH-60A helicopter at 20 ft.; network configuration 15-25-25-1; 

NCD (Pruned) learning rule. 
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Figure 18. Results for the UH-60A helicopter at 20 ft.; network configuration 15-25-25-1; 

NCD (Pruned) learning rule. 

(L 15-18-25-1 BPNN NCD 

Table 14 and Figure 19 present the results for this configuration. This 
setup is implemented as an alternative to the pruned network to reduce the network 
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training time and complexity. The RMS error is about 0.0384. The maximum error is 
2.27 knots at 35 kn ots. We note the error SD is larger at this speed. The overall error SD 


is 0.8853 knots, which is slightly higher than the NCD 15-25-25-1 network. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.8853 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.7156 

0.2597 

5.1943% 

1.7481 

15 

15.0498 

0.6162 

4.1081% 

1.3976 

25 

24.1842 

0.9682 

3.8726% 

1.7528 

35 

34.6029 

1.5548 

4.4424% 

2.2788 

45 

45.5511 

1.2541 

2.7870% 

2.2318 


Table 14. Results for the UH-60A helicopter at 20 ft.; network configuration 15-18-25-1; 

NCD learning rule. Results of UH-60A for 15-18-25-1 NCD. 
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Figure 19. Results for the UH-60A helicopter at 20 ft.; network configuration 15-18-25-1; 

NCD learning rule. 
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5, Baseline Data Set Analysis 

The baseline data set was formed combining both 20 ft. and 85 ft. data sets. First, 
the network was trained with this data set, and the network performance was measured 
using the baseline test set. Then, using single condition test data sets separately, the 
network performance was evaluated for IGE and OGE conditions. As a result, the 
network-training time increased, due to the increase in the number of data points in the 
baseline data set. 

a. 15-18-25-1BPNNNCD 

Results for this setup are shown in Table 15 and in Eigure 20. They show 
that the RMS error is 0.0755. The maximum error is about 6.3 knots when the speed is 35 
knots. The error SD of 1.4224 knots is higher than that observed with other baseline 
networks and the single condition data networks. The error percentage for all speeds at 1 
G is about 6 %. The error SD is higher at fast speed, especially when the sideslip angle is 
-60 and 60 degrees. 
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Eigure 20. Results for the EIH-60A helicopter with baseline data; network configuration 15- 


18-25-1; NCD learning rule. 
42 







Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.4224 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.9564 

0.3565 

7.1292% 

2.8415 

15 

16.4937 

0.9016 

6.0108% 

4.2300 

25 

26.0262 

1.5856 

6.3426% 

5.1133 

35 

38.8061 

2.3985 

6.8529% 

6.2953 

45 

45.4173 

2.1685 

4.8188% 

4.1447 


Table 15. Results for the UH-60A helicopter with baseline data; network configuration 15- 

18-25-1; NCD learning rule. 

b. 15-25-25-1BPNNNCD 

Results are depicted in Figure 21 and Table 16. They show that the RMS 
error is 0.0762, and the maximum error is 6.2424 knots at 45 kn ots. The airspeed error 
SD is larger when the speed is 25 knots. The largest error percentage is equal to 9 % (at 5 
knots). The overall network error SD is 1.5664 knots. 
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Figure 21. Results for the UF1-60A helicopter with baseline data; network configuration 15- 

25-25-1; NCD learning rule. 
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Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.5664 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.7732 

0.4005 

8.0108% 

1.8981 

15 

13.4327 

1.0494 

6.5960% 

3.1877 

25 

26.4841 

2.2515 

9.0058% 

6.1337 

35 

38.4891 

1.9183 

5.4810% 

5.1834 

45 

43.3588 

1.9560 

4.3466% 

6.2424 


Table 16. Results for the UH-60A helicopter with baseline data; network configuration 15- 

25-25-1; NCD learning rule. 

c. 15-25-25-1 BPNNExt. DBD 

Table 17 and in Figure 22 present the findings for this configuration. 
Results show that the RMS error is 0.0501, which is better than that of the NCD scheme. 
The maximum error is about 4.5 knots (at 35 knots). We note the predicted speed mean 
values are close to the target speeds. The error percentage at the speed of 5 knots is about 
10 %. The NN prediction is quite good at the speed of 45 knots. The overall network 
error SD is 0.7320, which is the best obtained with the baseline data set. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7320 

Mean of 

Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.2079 

0.5382 

10.7647% 

2.8105 

15 

16.2084 

0.6293 

4.1953% 

2.5304 

25 

25.0000 

0.7696 

3.0782% 

2.3555 

35 

37.3300 

1.0832 

3.0949% 

4.5313 

45 

45.9893 

0.8478 

1.8839% 

1.9112 


Table 17. Results for the UH-60A helicopter with baseline data; network configuration 15- 

25-25-1; Ext. DBD learning rule. 
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Figure 22. Results for the UFI-60A helicopter with baseline data; network configuration 15- 

25-25-1; Ext. DBD learning rule. 


d. 15-25-25-1 BPNNNCD (Pruned) 

Table 18 and Figure 23 present the findings for this configuration. Results 
show that one PE was pruned in the first hidden layer. The resulting RMS error is 0.1213, 
and the maximum error at 45 knots is 8.6414 knots. The largest error SD and error 
percentage (11.5 %) are obtained for a 25 knots speed. The overall network error SD is 
2.0830 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 2.0830 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.9069 

0.3420 

6.8404% 

2.8877 

15 

9.8217 

1.1116 

7.4108% 

7.2736 

25 

23.5765 

2.8864 

11.5454% 

6.5795 

35 

38.3384 

2.8346 

6.6262% 

7.1468 

45 

43.8399 

2.9818 

1.8839% 

8.6414 


Table 18. 


Results for the E1F1-60A helicopter with baseline data; network configuration 15- 


25-25-1; NCD (Pruned) learning rule. 


45 

























Predicted Speed vs Grossweight 


Predicted Speed vs Actual Speed 



1.6 1.8 2 2.2 2.4 

Gross Weight (lb* 10000) 


Slow Speed vs Angle 


30 

£25 

■o 

Of 

4f 

O. 

</> 

■S 15 

U 

s 

CL 

5 

0 


.. 

eeoeeooooooo 

0 60 120 180 240 300 360 



Anual Speed (l<^) 

Fast Speed vs Angle 



Side Slip A^gle (Deg) 


Side Slip A^gle (Deg) 


15-25-25-1 BPNN NCD (Prune) UH-60A RMS = 0.1213 ALT = 20 ft and 85 ft (Sea Level) 


Figure 23. Results for the UH-60A helicopter with baseline data; network configuration 15- 

25-25-1; NCD (Pruned) learning rule. 
e. Baseline Data Set IGE Analysis 

(1) 15-25-25-1 BPNN NCD. Results are shown in Table 19 

and Figure 24. They show the RMS error is 0.0759. The maximum error is 6.2 knots, 
which is almost the same for speeds 25 knots and higher. The network performance at 
high speeds is better than that of at low speeds since the percent error at 35 and 45 knots 
is about 5 %, whereas at lower speeds it is twice that number. The overall error SD is 
1.5516 knots. 


Actual 

Airspeed (kts) 

NEURAL N] 
Tota 

ETWORK RESULTS 
[SD = 1.5516 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.7699 

0.5063 

10.1268% 

1.8981 

15 

13.7003 

1.1242 

7.4947% 

2.7778 

25 

26.5671 

2.1650 

8.6599% 

6.1337 

35 

38.3892 

1.7970 

5.1342% 

5.1715 

45 

42.8593 

2.0669 

4.5931% 

6.2424 


Table 19. Results for the UH-60A helicopter at 20 ft with baseline data; network 


configuration 15-25-25-1; NCD learning rule. 
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Figure 24. Results for the UH-60A helieopter at 20 ft with baseline data; network 
configuration 15-25-25-1; NCD learning rule. 

(2) 15-25-25-1 BPNN Ext, DBD, Results are shown in Table 

20 and Figure 25. They show the RMS error is 0.0499. The maximum error of 4.5 knots 


is obtained at a 35 knots speed. Note that the network performance is quite good when 
compared with the NCD scheme. The overall error SD is 0.6505 knots, which is 


significantly less than it is in the NCD scheme. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.6505 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.5127 

0.4383 

8.7656% 

2.2475 

15 

16.4063 

0.6193 

4.1287% 

2.5304 

25 

25.1852 

0.7545 

3.0179% 

2.3505 

35 

37.6802 

0.9720 

2.7772% 

4.5313 

45 

46.0822 

0.4932 

1.0961% 

1.8867 


Table 20. Results for the UH-60A helicopter at 20 ft with baseline data; network 

configuration 15-25-25-1; Ext. DBD learning rule. 
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Figure 25. Results for the UH-60A helicopter at 20 ft with baseline data; network 
configuration 15-25-25-1; Ext. DBD learning rule. 


(3) 15-25-25-1 BPNN NCD (Pruned). Results are shown in 

Table 21 and Figure 26. They show that the RMS error is 0.0901. The maximum speed 
error of 7.9 knots is obtained for a 45 knots speed. We note that the network performance 


is not good when the pruning facility is enabled. The overall error SD is 1.7173 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.7173 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

4.3654 

0.2933 

5.8653% 

1.1256 

15 

11.5338 

1.1787 

7.8578% 

4.9166 

25 

23.9066 

2.2817 

9.1268% 

4.0707 

35 

37.3058 

2.1850 

6.2429% 

4.4207 

45 

42.1901 

2.5806 

5.7346% 

7.9128 


Table 21. Results for the UH-60A helicopter at 20 ft with baseline data; network 
configuration 15-25-25-1; NCD (Pruned) learning rule. 
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Figure 26. Results for the UH-60A helicopter at 20 ft with baseline data; network 
configuration 15-25-25-1; NCD (Pruned) learning rule. 
f. Baseline Data Set OGE Analysis 

(1) 15-25-25-1 BPNN NCD, Findings for this configuration 

are presented in Table 22 and Figure 27. Results show that the RMS error is 0.0766. The 
maximum error is 5.86 kn ots. The network performance at a speed of 25 knots is not 
good, as the error SD is 2.3 knots at that speed. The error percentage is about 5 % for 


speeds other than 25 knots but it increases up to 9 % for 25 knots. The overall error SD is 
1.5516 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.5516 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.7779 

0.2569 

5.1390% 

1.6144 

15 

13.1631 

0.8966 

5.9771% 

3.1890 

25 

26.4013 

2.3416 

9.3664% 

5.8619 

35 

38.5880 

2.0476 

5.8502% 

5.1844 

45 

43.8585 

1.7192 

3.8204% 

4.3961 


Table 22. Results for the UH-60A helicopter at 85 ft with baseline data; network 


configuration 15-25-25-1; NCD learning rule. 
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Figure 27. Results for the UH-60A helicopter at 85 ft with baseline data; network 
configuration 15-25-25-1; NCD learning rule. 

(2) 15-25-25-1 BPNN Ext. DBD. Table 23 and in Figure 28 

present the results for this configuration. Results show that the RMS error is 0.0482. The 
maximum error observed at 35 knots is equal to 3.8 kn ots. Note that this error is smaller 
than that observed in the NCD scheme, and the network performance is quite good when 
compared with the NCD scheme for the OGE baseline setup. At low speeds, the error SD 
is about 0.5 knots. The network performance is quite good for all speeds except for 5 
knots where the error percentage increases to 9 %. The overall error SD is 0.6850 knots. 
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Figure 28. Results for the UH-60A helicopter at 85 ft with baseline data; network 
configuration 15-25-25-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.6850 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.9227 

0.4746 

9.5519% 

2.8377 

15 

15.8606 

0.5962 

3.9746% 

1.9213 

25 

24.8085 

0.6801 

2.7205% 

1.5601 

35 

36.8469 

1.0210 

2.9171% 

3.8006 

45 

45.9612 

0.9043 

2.0096% 

1.9007 


Table 23. Results for the UH-60A helicopter at 85 ft with baseline data; network 

configuration 15-25-25-1; NCD learning rule. 

(3) 15-25-25-1 BPNN NCD (Pruned). Table 24 and in 

Figure 29 present the results for this configuration. Results show that the RMS error is 
0.1087. The maximum error observed for a 35 knots speed is equal to 7.14 knots. Note 
that the network performance degrades when the pruning facility is enabled. The overall 
error SD is 1.8642 knots. 
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Figure 29. Results for the UH-60A helicopter at 85 ft with baseline data; network 
configuration 15-25-25-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.8642 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.7243 

0.1567 

3.1338% 

1.6318 

15 

10.3127 

0.5639 

3.7592% 

5.7803 

25 

24.6236 

2.3566 

9.4263% 

4.6493 

35 

39.1423 

3.2081 

9.1660% 

7.1416 

45 

45.3547 

2.8266 

6.2814% 

3.5760 


Table 24. Results for the UH-60A helicopter at 85 ft with baseline data; network 
configuration 15-25-25-1; NCD (Pruned) learning rule. 

6, Simplifying The Data Set Using Eigenvalues and Eigenvectors 

Preprocessing the data to obtain a simplified NN is often useful. Goff developed a 
NN helicopter airspeed prediction study similar to the Haas and McCool study [Ref. 20]. 
Goff analyzed the each input contribution to the network performance and determined 
that 11 of the 16 inputs played an important role. Results obtained by preprocessing the 
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data and simplifying the network inputs were quite suceessful. His study on the Lynx 
MK9 showed that airspeed in the low speed environment eould be predieted with an error 
of +/-3.1 knots at 1 a [Ref 20]. 

Another way to deerease the NN input dimension and simplify the network 
strueture by projeeting the high-dimensional data onto a lower dimensional input spaee 
using prineipal eomponent analysis (PCA). PCA is applied by eonsidering the input data 
set as a matrix, A, and translating this matrix to a diagonal or upper triangular form to 
eompute its eigenvalues and eigenveetors. Eigenveetors are the normal modes of the 
system and they aet independently. The beauty of eigenvalues and eigenveetors is that 
they eapture the eharaeteristies and behavior of the whole system. The key equation for 
eigenvalues and eigenveetors is: 

Ax=Xx, (4.1) 

where k is an eigenvalue of the matrix A and the veetor x is the assoeiated eigenveetor. 
Most veetors x will not satisfy this equation. A typieal x ehanges direetion when 
multiplied by A, so that Ax is not a multiple ofx. Thus, only eertain speeial numbers are 
eigenvalues and only eertain speeial veetors are eigenveetors [Ref 18]. 

After obtaining the eigenvalues and eigenveetors of the single eondition data set, 
dominating eigenvalues were seleeted. It was observed that six of the eigenvalues were 
the dominating eigenvalues, leading to a input data sub-matrix with a dimension of 
[14x6]. This sub-matrix was multiplied by the whole data set again. The purpose of this 
proeedure was to reduee the number of inputs to 6 while still eapturing the properties 
from all inputs. Several network arehiteetures were implemented using this simplified 
data set. Network eonfigurations with the best performanees are deseribed in the 
following seetions. Note that none of these performed as well as the best network trained 
with the full data sets. 

a. 6-18-18-1 BPNNExt DBD 

Results of this setup are shown in Figure 30 and in Table 25. They show 
that the RMS error is 0.0565. Although the maximum error of 5.4 knots is at 45 knots, the 
11% error pereentage is higher for a 5 knots speed. The overall error SD is 1.1692 knots. 
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The network performance at low speeds is not as satisfactory as it is at fast speeds since 
the percent error is quite large at low speeds. 
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Figure 30. Results for the UF1-60A helicopter at 85 ft with simplified data; network 


configuration 6-18-18-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.1692 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.9375 

0.5879 

11.7571% 

3.0269 

15 

15.3870 

0.7861 

5.2444% 

2.2.1139 

25 

23.8447 

1.2453 

4.9813% 

3.5517 

35 

35.4789 

1.4587 

4.1677% 

2.5422 

45 

43.5555 

2.2160 

4.9244% 

5.3797 


Table 25. Results for the UH-60A helicopter at 85 ft with simplified data; network 


configuration 6-18-18-1; Ext. DBD learning rule. 

b. 6-18-18-1 BPNNNCD (Pruned) 

Results are depicted in Eigure 31 and Table 26. One PE in the first layer 

was disabled by the NN. Results show that the RMS error is 0.0648. The maximum error 

is over 5.0 knots at speeds of 25 and 35 knots. The error SD is also larger at those speeds. 

We note that, the maximum error percentage is 7.95 % at a speed of 5 knots, and the 
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overall error SD is 1.1942 knots. In conclusion, results showed that the network 
performance is not as good as that obtained when using the NCD learning rule without 
the pruning facility. 
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Figure 31. Results for the UH-60A helicopter at 85 ft with simplified data; network 

configuration 6-18-18-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

:twork results 

[SD = 1.1942 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.5084 

0.3977 

7.9568% 

1.9809 

15 

14.7297 

0.5445 

3.6297% 

1.1050 

25 

26.7339 

1.8057 

7.2229% 

5.4495 

35 

38.2674 

1.8897 

5.3992% 

5.3390 

45 

43.6582 

0.9459 

2.1020% 

2.6491 


Table 26. Results for the UH-60A helicopter at 85 ft with simplified data; network 
configuration 6-18-18-1; NCD (Pruned) learning rule. 
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c. 6-18-18-1BPNNNCD 

Figure 32 and Table 27 present the results for this eonfiguration. Results 
show that the RMS error equal to 0.06561 is the best among all the configurations 
investigated with the simplified data. The maximum error is about 4 knots when the 
helicopter is moving at a speed of 35 knots with a 30*’ sideslip angle to the left. The error 
SD is above 1.5 knots at that speed. The error percentage for each speed is about 5 %, 
except for a 45 knots speed, and the overall error SD is equal to 1.0032 knots. 
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Figure 32. Results for the UH-60A helicopter at 85 ft with simplified data; network 


configuration 6-18-18-1; NCD learning rule 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.0032 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.7558 

0.2888 

5.5765% 

2.7471 

15 

15.3714 

0.6022 

4.0149% 

1.9253 

25 

25.0220 

1.3322 

5.3287% 

2.7581 

35 

37.1657 

1.5913 

4.5467% 

3.8912 

45 

43.5339 

1.2950 

2.8777% 

3.1240 


Table 27. Results for the UF1-60A helicopter at 85 ft with simplified data; network 


configuration 6-18-18-1; NCD learning rule. 
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Results showed that the network with simplified data performs worse than 
the network with the pruning faeility enabled. Therefore, the simplified data was not used 
for OH-6A analyses. 

B, ANALYSIS OF THE OH-6A MODEL 

1. Out of Ground Effect Analysis at Sea Level 

This seetion presents results of the OGE analysis for the OH-6A helieopter with 
the simulator altimeter set at 100 feet AGE. The network was trained using the single 
eondition data set obtained for this altitude similarly to the EIH-60A analysis. Results for 
the OGE analysis using the eombined data set are shown in the OH-6A Baseline Analysis 
Section. Note that the engine torque input parameter for the NN model was removed due 
to the limitation of the engine model parameters for the OH-6A helicopter in the 
EEIGHTEAB simulator. Therefore, a 14-25-25-1 setup was used for all OH-6 A model 
analyses. Einally, in light of the results obtained from the EIH-60A helicopter analysis, 
only the NCD and Ext. DBD learning rules were investigated. 
a. 14-25-25-1BPNNNCD 

Results for this setup are given in Eigure 33 and Table 28. They show that 
the RMS error is 0.0609. The 5 knots maximum speed error is observed at 35 knots, 
while the absolute maximum error is less than 3 knots at other speeds. The airspeed error 
SD is less than 1 knot except for 45 knots. The overall error SD is 0.7921 knots. 


Actual 

Airspeed (kts) 

NEURAE NETWORK RESUETS 

Total SD = 0.7921 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

4.2398 

0.4797 

9.5947% 

1.5363 

15 

12.7387 

0.7130 

4.7531% 

3.3062 

25 

25.6227 

0.8655 

3.4618% 

2.0925 

35 

38.6028 

0.8703 

2.4865% 

5.2048 

45 

44.8812 

1.2253 

2.7229% 

1.9604 


Table 28. Results for the OH-6A helicopter at 100 ft (SE); network configuration 

14-25-25-1; NCD learning rule. 
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Figure 33. Results for the OH-6A helieopter at 100 ft (SL); network configuration 

14-25-25-1; NCD learning rule. 

b. 14-25-25-1 BPNNExt DBD 

Results are shown in Figure 34 and Table 29. They show that the RMS 
error is 0.0735. Note that the absolute maximum error is about 5 knots at the speed of 35 
knots, but the overall error SD of 1.3207 knots is larger than that observed for the NCD 


setup. For all speeds, except for 45 knots, the error percentage is 5 % and over. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.3207 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.6946 

0.4724 

9.4478% 

2.3008 

15 

13.3800 

1.2846 

8.5642% 

3.8498 

25 

25.3109 

1.2145 

4.8578% 

3.2181 

35 

36.3450 

2.6079 

7.7543% 

5.0472 

45 

43.8989 

1.1780 

2.6179% 

3.2799 


Table 29. Results for the OH-6A helicopter at 100 ft (SL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
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Figure 34. Results for the OH-6A helicopter at 100 ft (SL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
c. 14-25-25-1 BPNNNCD (Pruned) 

Figure 35 and Table 30 present the findings for this configuration. Results 
show that the RMS error is 0.0669. The maximum speed error is about 4.6 knots at a 
speed of 35 knots. The error SD is less than 0.88 knots except for a 45 knots speed. The 
network performance is quite good except for 5 knots as the error percentage is equal to 9 


%, which is very large at that speed. The overall error SD is 0.759 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.7590 

Mean of 

Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.9667 

0.4976 

9.9521% 

1.9224 

15 

12.0952 

0.6239 

4.1595% 

4.0568 

25 

25.1773 

0.8877 

3.5509% 

1.7791 

35 

38.4750 

0.6479 

1.8510% 

4.6989 

45 

44.9067 

1.2291 

2.7914% 

1.9686 


Table 30. Results for the OH-6A helicopter at 100 ft (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 
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Predicted Speed vs Grosswreight 


Predicted Speed vs Actual Speed 
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Figure 35. Results for the OH-6A helicopter at 100 ft (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 

2, In-Ground Effect Analysis at Sea Level 

The OH-6 A helicopter in-ground effect analysis was performed at 12 ft AGL. 
Same network architectures as those considered with out-of-ground analyses were 
considered. Results are presented below. 

a. 14-25-25-1 BPNN NCD 

Results are given in Table 31 and Figure 36. They show that the RMS 
error for this architecture is 0.0592. The maximum error is 5 knots for a speed equal to 35 
knots. The percentage of error at 1 g is about 3% for all speeds except for 5 kn ots. The 
speed error SD gets larger at speeds above 25 kn ots. The overall error SD is 0.8465 
knots. 
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Figure 36. Results for the 0H-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; NCD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.8465 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.4013 

0.4882 

9.7633% 

1.4300 

15 

13.0112 

0.3802 

2.5348% 

2.7769 

25 

25.9480 

1.1424 

4.5695% 

2.7853 

35 

38.5135 

0.9581 

2.7373% 

5.1048 

45 

45.3026 

1.2993 

2.8872% 

2.1848 


Table 31. Results for the OF1-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; NCD learning. 


b. 14-25-25-1 BPNN Ext DBD 

Table 32 and Figure 37 present the findings for this configuration. Results 

show that the RMS error is 0.062. The maximum error is equal to 6 knots (at 35 knots). 

The maximum error percentage is 12 % and occurs at 5 knots. The network performance 

is quite good at 45 knots with an error SD of 0.8 knots. We note that the overall error SD 

of 0.9678 is slightly higher than that of obtained with the NCD scheme. 
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Actual 

Airspeed (kts) 

NEURAL N] 
Tota 

ETWORK RESULTS 

SD = 0.9678 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.9451 

0.6122 

12.2449% 

1.1009 

15 

15.4143 

0.6925 

4.6168% 

1.6649 

25 

27.1087 

1.1644 

4.6577% 

3.7000 

35 

38.7739 

1.6571 

4.7345% 

5.9939 

45 

45.7621 

0.8066 

1.7952% 

2.0848 


Table 32. Results for the OH-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
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Figure 37. Results for the OH-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; Ext. DBD learning rule, 
c. 14-25-25-1 BPNNNCD (Pruned) 

Figure 38 and Table 33 present the results for this setup. Results show that 
the RMS error is 0.0674 and the maximum error of 4.4 knots occurs at 35 knots. Note 
that the maximum error is about 2 knots for all other speeds. The error SD is the largest at 
25 and 45 kn ots. The overall network error SD is 0.8803 knots. 
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Actual 
Airspeed (kts) 

NEURAL N 
Tota 

ETWORK RESULTS 

SD = 0.8803 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

4.1006 

0.6088 

12.1766% 

1.9135 

15 

12.2608 

0.3937 

2.6246% 

3.4638 

25 

25.5791 

1.1830 

4.7320% 

2.5694 

35 

38.5929 

0.5935 

1.6956% 

4.4157 

45 

45.1430 

1.4917 

3.3148% 

2.0239 


Table 33. Results for the OH-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 
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Figure 38. Results for the OH-6A helicopter at 12 ft (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 

3, OH-6A Baseline Data Analysis at Sea Level 

The baseline data set was formed by combining both 12 ft. and 100 ft. data sets. 
First the network was trained with this combined data set, and the network performance 
was measured with the baseline test set. Then, using single condition test data sets 
separately, the network performance was evaluated for IGF and OGE conditions. 
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a. 14-25-25-1BPNNNCD 

Results for this configuration are shown in Figure 39 and Table 34. They 
show that the RMS error is 0.0622. The absolute maximum speed error occurs at 35 knots 
and is equal to 4.82 knots. The airspeed error SD at 1 a is more than 1 knot at the speed 
of 25 knots and over. The overall error SD is 1.2 knots. 
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Figure 39. Results for the OH-6A helicopter with baseline data (SL); network configuration 


14-25-25-1; NCD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESUETS 

Total SD = 1.1926 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.6199 

0.4348 

8.6964% 

2.3155 

15 

13.2606 

0.4590 

3.0599% 

2.9471 

25 

25.8467 

1.6476 

6.5904% 

4.3737 

35 

37.8501 

1.1617 

3.3190% 

4.8224 

45 

44.8255 

2.0805 

4.6234% 

2.8179 


Table 34. Results for the OH-6A helicopter with baseline data (SL); network configuration 

14-25-25-1; NCD learning rule. 
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b. 14-25-25-1 BPNNExt. DBD 

Table 35 and Figure 40 present the results for this setup. Results show that 


the RMS error is 0.0687. The maximum speed error observed is 7.5 knots for 35 kn ots. 
At most speeds, the error SD at 1 a is more than 1.2 knots. Hence, the error percentages 
are also larger when compared with other network performances. The overall error SD is 
1.4071 knots. 
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Figure 40. Results for the OH-6A helicopter with baseline data (SL); network configuration 


14-25-25-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL N] 
Tota 

ETWORK RESULTS 
[SD = 1.4071 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.9860 

0.7204 

14.4071% 

3.5396 

15 

14.7369 

1.3806 

9.2038% 

4.0839 

25 

26.0708 

1.2023 

4.8092% 

4.9530 

35 

38.3260 

2.6810 

7.6601% 

7.4491 

45 

45.6813 

1.3554 

3.0120% 

2.9237 


Table 35. Results for the OH-6A helicopter with baseline data (SL); network configuration 

14-25-25-1; Ext. DBD learning rule. 


65 









































c. 14-25-25-1 BPNNNCD (Pruned) 

Findings for this configuration are presented in Table 36 and Figure 41. 
The pruning disabled 3 PEs in the first hidden layer. The RMS error is 0.0611 and the 
maximum speed error is 4.98 knots at 25 knots. The error SD is greater than 1.2 knots for 
the speeds 25 knots and higher. The overall error SD is 1.1689 knots. 
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Figure 41. Results for the OF1-6A helicopter with baseline data (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.1689 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.9731 

0.3442 

6.8035% 

1.7741 

15 

13.4709 

0.3862 

2.5744% 

2.3807 

25 

26.1026 

1.6586 

6.6343% 

5.0413 

35 

38.1568 

1.2776 

3.6504% 

4.9842 

45 

44.9858 

1.9399 

4.3109% 

3.1513 


Table 36. Results for the OF1-6A helicopter with baseline data (SL); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 
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d. 0H-6A SL Baseline Data IGE Analysis 

(1) 14-25-25-1 NCD IGE, Table 37 and Figure 42 present 

detailed results for this configuration. Results show that the RMS error is 0.0593. The 
maximum speed error is 4.5 knots observed for a 35 knots speed. The largest error SD is 
2.14 knots at 45 knots. The overall error SD is 1.1501 knots. 
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Figure 42. Results for the OH-6A helicopter at 12 ft with baseline data (SL); network 

configuration 14-25-25-1; NCD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESUETS 

Total SD = 1.1501 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.7499 

0.4647 

9.2946% 

2.2505 

15 

13.4738 

0.3264 

2.1762% 

2.1856 

25 

26.0936 

1.5397 

6.1587% 

4.2902 

35 

37.6622 

1.1159 

3.1883% 

4.4943 

45 

44.6424 

2.1425 

4.7612% 

2.8106 


Table 37. Results for the OH-6A helicopter at 12 ft with baseline data (SL); network 

configuration 14-25-25-1; NCD learning rule. 


67 
































(2) 14-25-25-1 NCD (Pruned) IGE. Results are given in 

Table 38 and Figure 43. They show that the RMS error is 0.0592. The maximum speed 
error is 4.5 knots observed for a 35 knots speed. The overall error SD is 1.155 knots. We 
note that no significant improvement was obtained when compared with the un-pruned 
NCD scheme. 
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Figure 43. Results for the OH-6A helicopter at 12 ft with baseline data (SL); network 


configuration 14-25-25-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESUETS 

Total SD = 1.1549 

Mean of 

Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

At 1 G 

Abs. Maximum 
Error (kts) 

5 

3.8734 

0.4085 

8.1709% 

2.0251 

15 

13.5302 

0.3172 

2.1143% 

2.0447 

25 

26.1493 

1.5478 

6.1914% 

4.1953 

35 

37.7847 

1.0744 

3.0696% 

4.5171 

45 

44.7929 

2.2042 

4.8981% 

2.8103 


Table 38. Results for the OH-6A helicopter at 12 ft with baseline data (SL); network 
configuration 14-25-25-1; NCD (Pruned) learning rule. 
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(3) 14-25-25-1 Ext, DBD IGE, Results are illustrated in 

Figure 44 and Table 39. They show that the RMS error is 0.0686. The maximum speed 
error is 7.4491 knots observed for a 35 knots speed. We note that neither the prediction of 
low speed nor the prediction of fast speed is as good as that obtained with the NCD setup. 
The overall error SD is 1.3293 knots. 
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Figure 44. Results for the OF1-6A helicopter at 12 ft with baseline data (SL); network 


configuration 14-25-25-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.3293 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.4696 

0.5957 

11.9134% 

2.7024 

15 

15.3847 

1.1166 

7.4493% 

2.5874 

25 

26.5757 

1.1507 

4.6029% 

4.9330 

35 

38.6624 

2.7151 

7.7574% 

7.4491 

45 

45.8415 

1.4335 

3.1855% 

2.9237 


Table 39. Results for the OF1-6A helicopter at 12 ft with baseline data (SL); network 

configuration 14-25-25-1; Ext. DBD learning rule. 
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e. 0H-6A SL Baseline Data OGE Analysis 

(1) 14-25-25-1 Ext, DBD OGE, Results for this network are 

depicted in Figure 45 and Table 40. They show that the RMS error is 0.0687. The 
maximum speed error is 6.6 knots observed for a 35 knots speed. The error SD is over 1 
knot almost at all speeds. The overall error SD is 1.3015 knots. 
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Figure 45. Results for the OH-6A helicopter at 100 ft with baseline data (SL); network 


configuration 14-25-25-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD= 1.3015 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.5071 

0.4741 

9.4826% 

3.5396 

15 

14.0892 

1.3168 

8.7787% 

4.0839 

25 

25.5659 

1.0316 

4.1264% 

3.7463 

35 

37.9896 

2.6280 

7.5087% 

6.6395 

45 

45.5210 

1.2654 

2.8119% 

2.2839 


Table 40. Results for the OH-6A helicopter at 100 ft with baseline data (SL); network 

configuration 14-25-25-1; Ext. DBD learning rule. 
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(2) 14-25-25-1 BPNN NCD OGE. Figure 46 and Table 41 

present the results for this setup. Results show that the RMS error is 0.0671. The 
maximum speed error is 4.8 knots observed for a 35 knots speed. The error SD is above 1 
knot for speeds 25 knots and higher. The overall error SD is 1.2022 knots. 
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Figure 46. Results for the OH-6A helicopter at 100 ft with baseline data (SL); network 

configuration 14-25-25-1; NCD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 1.2022 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.4898 

0.3599 

7.1982% 

2.3155 

15 

13.0474 

0.4741 

3.1618% 

2.9471 

25 

25.5998 

1.7194 

6.8786% 

4.3737 

35 

38.038 

1.1861 

3.3887% 

4.8224 

45 

45.0086 

2.0195 

4.4874% 

2.8176 


Table 41. Results for the OH-6A helicopter at 100 ft with baseline data (SL); network 

configuration 14-25-25-1; NCD learning rule. 
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(3) 14-25-25-1 BPNN NCD (Pruned) OGE. Results are 

depicted in Figure 47 and in Table 42. Results show that the RMS error is 0.0650. The 
maximum speed error is about 5 knots at 35 knots. The maximum error is about 2 knots 
for all other speeds. One PE in the first hidden layer was disabled. The overall error SD is 
1.2 knots. Note that performance of this setup is better than that of the Ext. DBD 
configuration. 
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Eigure 47. Results for the OH-6A helicopter at 100 ft with baseline data (SE); network 


configuration 14-25-25-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAE NETWORK RESUETS 

Total SD= 1.2096 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.8439 

0.3031 

6.0621% 

1.8782 

15 

13.3969 

0.5177 

3.4513% 

2.5735 

25 

25.9133 

1.7329 

6.9317% 

4.5386 

35 

38.2975 

1.184 

3.3828% 

5.1005 

45 

45.2226 

2.0303 

4.5118% 

2.8686 


Table 42. Results for the OH-6A helicopter at 100 ft with baseline data (SE); network 


configuration 14-25-25-1; NCD (Pruned) learning rule. 
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4, OH-6A Out of Ground Effect Analysis at High Altitude 

The simulator was run at a pressure altitude of 6000 feet for high altitude analysis. 
At that pressure altitude, the altitude AGL was set to 100 ft. for the OGE condition, 
which is the same AGL altitude, that the sea level OGE analysis was performed at. 

a. 14-25-25-1BPNNNCD 

Results are presented in Eigure 48 and Table 43. They show that the RMS 
error for this setup is 0.0485. The maximum speed error, which is 2.47 knots, and the 
maximum error percentage equal to 9.4 % are both observed for a 5 knots speed. The 
maximum airspeed error SD of 1.1356 knots occurs at 35 knots. The overall error SD is 
0.6637 knots. 
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Eigure 48. Results for the OH-6A helicopter at 100 ft (HE); network configuration 

14-25-25-1; NCD learning rule. 
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Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.6637 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.5101 

0.4738 

9.4794% 

2.4762 

15 

16.209 

0.6861 

4.5744% 

2.2914 

25 

23.272 

0.5601 

2.2405% 

2.3443 

35 

35.288 

1.1356 

3.2444% 

1.8547 

45 

45.985 

0.6179 

1.3732% 

1.8831 


Table 43. Results for the OH-6A helicopter at 100 ft (HL); network configuration 

14-25-25-1; NCD learning rule. 


b. 14-25-25-1 BPNNExt DBD 

Figure 49 and Table 44 present the findings for this configuration. Results 
show that the RMS error is 0.0958. This is one of the largest errors observed for the OH- 
6A analysis. The maximum error is 6.36 knots, which is more than twice the error of the 
NCD setup. The error percentage at 5 knots is better than that of the NCD configuration. 
The overall error SD is 1.022 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.022 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.1676 

0.2838 

5.6769% 

2.3937 

15 

19.875 

0.7859 

5.2395% 

6.3677 

25 

23.339 

1.2915 

5.1659% 

5.3069 

35 

33.995 

1.9872 

5.6777% 

3.4085 

45 

45.699 

0.3127 

0.6917% 

1.1757 


Table 44. Results for the OH-6A helicopter at 100 ft (HL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
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Figure 49. Results for the 0H-6A helieopter at 100 ft (HL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
c. 14-25-25-1 BPNNNCD (Pruned) 

Figure 50 and Table 45 present the findings for this configuration. The 
pruning disabled 1 PE of the first hidden layer. Results show that the RMS error is 
0.0475. The maximum speed error is 2.85 knots observed for 5 knots speed. The results 
for fast speed predictions are better than the those at low speeds. The overall error SD is 
0.6401 knots. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.6401 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.8147 

0.3546 

7.0932% 

2.8415 

15 

14.812 

0.7660 

5.1070% 

1.4709 

25 

24.0 

0.7102 

2.8410% 

2.2145 

35 

36.208 

0.7016 

2.0047% 

2.0621 

45 

46.040 

0.6219 

1.3821% 

2.0527 


Table 45. Results for the OH-6A helicopter at 100 ft (HE); network configuration 

14-25-25-1; NDC (Pruned) learning rule. 
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Figure 50. Results for the OH-6A helicopter at 100 ft (HL); network configuration 

14-25-25-1; NDC (Pruned) learning rule. 


5, OH-6A In-Ground Effect Analysis at High Altitude 

In-ground effect analyses of the OH-6A at high altitude were implemented by 
setting the simulator pressure altitude to 6000 feet. At that pressure altitude, the altitude 
AGL was set to 12 ft. for the IGE condition, which is the same AGL altitude at which sea 
level IGE analyses were performed. 

a. 14-25-25-1 BPNN NCD 

Eigure 51 and Table 46 present the findings for this configuration. Results 
show that the RMS error is 0.0516, and the maximum speed error is 2.596 knots for a 5 
knots speed. The maximum error percentage equal to 9.12 % is also observed at that 
speed. The maximum airspeed error SD of 1.20 knots occurs at 35 kn ots. The overall 
error SD is 0.6855 knots. 
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Figure 51. Results for the OH-6A helicopter at 12 ft (HL); network configuration 

14-25-25-1; NDC learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESUETS 

Total SD = 0.6855 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.4367 

0.4811 

9.1210% 

2.5962 

15 

16.3790 

0.6755 

4.503% 

2.4700 

25 

23.1736 

0.5637 

2.2547% 

2.5948 

35 

35.3076 

1.2062 

3.4462% 

2.0407 

45 

46.0330 

0.7075 

1.5722% 

2.0734 


Table 46. Results for the OH-6A helicopter at 12 ft (HL); network configuration 

14-25-25-1; NDC learning rule. 


b. 14-25-25-1 BPNN Ext DBD 

Results for this configuration are presented in Figure 52 and Table 47. 
They show that the RMS error is 0.0865, which is larger than that of the NCD setup. 
Although the maximum error is 5.84 knots for a 15 knots speed, the maximum error 
percentage equal to 5.53 % is significantly smaller than that of the NCD scheme. The 
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maximum airspeed error SD of 1.71 knots oceurs at 35 knots. Note that the airspeed 
estimation is poor at low speeds and it is best at high speeds when the gross weight of the 
helieopter is less than 2100 lb. The overall error SD is 0.9556 knots. 
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Figure 52. Results for the OH-6A helieopter at 12 ft (HL); network configuration 

14-25-25-1; Ext. DBD learning rule. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.9556 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.5658 

0.2767 

5.5336% 

1.9728 

15 

19.5361 

0.8258 

5.5054% 

5.8437 

25 

23.7345 

1.1813 

4.7250% 

3.7400 

35 

34.4205 

1.7182 

4.9093% 

2.8189 

45 

45.5856 

0.6026 

1.3391% 

1.6418 


Table 47. Results for the OF1-6A helicopter at 12 ft (FIL); network configuration 

14-25-25-1; Ext. DBD learning rule. 
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c. 14-25-25-1 BPNNNCD (Pruned) 

Results are presented in Figure 53 and Table 48. One PE in the first 
hidden layer was disabled. Results show that the RMS error is 0.0443, which is the best 
result obtained for the analysis of this altitude. The maximum speed error is 2.58 knots 
(at 5 knots) and the maximum error percentage is 5.59 %. The maximum airspeed error 
SD equal to 0.5883 knots occurs at 15 knots speed. The overall error SD is 0.6401 knots. 
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Figure 53. Results for the OH-6A helicopter at 12 ft (HE); network configuration 


14-25-25-1; NCD (Pruned) learning rule. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.443 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

2.9434 

0.2797 

5.5945% 

2.5801 

15 

14.7437 

0.6863 

4.5756% 

1.4916 

25 

24.2145 

0.6529 

2.6118% 

2.0764 

35 

36.3517 

0.6982 

1.9949% 

2.1910 

45 

46.8740 

0.6187 

1.3749% 

1.8320 


Table 48. Results for the OF1-6A helicopter at 12 ft (FIE); network configuration 

14-25-25-1; NCD (Pruned) learning rule. 
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6, OH-6A Baseline Data Analysis at High Level 

A baseline training data set was formed by eombining both 12 ft. and 100 ft. set of 
high level altitude data. First, the network was trained with this combined data set, and 
then the network performance was measured using the baseline testing set. Finally, the 
network performance was evaluated for IGF and OGE conditions, using single condition 
test data sets separately. Results are given below. 

a. 14-25-25-1BPNNNCD 

Results are shown in Figure 54 and Table 49. They show that the RMS 
error is 0.0507, and the absolute maximum speed error is 3.15 knots (at 25 knots). The 
airspeed error SD at 1 a is less than 1 knot at all speeds. The overall error SD is 0.7139 
knots. 
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Figure 54. Results for the OH-6A helicopter with baseline data (HE); network configuration 

14-25-25-1; NCD learning rule. 
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Aetual 

Airspeed (kts) 

NEURAL NI 
Total 

iTWORK RESULTS 

1SD = 0.7139 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Pereent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.1503 

0.4366 

8.7318% 

2.7569 

15 

15.9246 

0.8179 

5.4525% 

2.1210 

25 

23.3005 

0.7533 

3.0132% 

3.1582 

35 

35.7437 

0.9455 

2.7014% 

2.5005 

45 

45.8256 

0.6244 

1.3876% 

2.1869 


Table 49. Results for the OH-6A helieopter with baseline data (HL); network eonfiguration 

14-25-25-1; NCD learning rule. 


b. 14-25-25-1 BPNNExt DBD 

Results are presented in Figure 55 and Table 50. They show that the RMS 
error is 0.0772, and the absolute maximum speed error is 5.39 knots, whieh oeeurs at 15 
knots speed. The maximum airspeed error SD at 1 a is 1.61 knots (at 25 knots). The 
predietion aeeuraey of low speeds is not as good as that of fast speeds. The overall error 
SD is 1.0209 knots. 


Aetual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.0209 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Pereent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.7873 

0.4398 

8.7951% 

2.0195 

15 

19.0707 

0.7415 

4.9436% 

5.3927 

25 

24.8423 

1.6157 

6.4628% 

3.7711 

35 

35.460 

0.9531 

2.7230% 

2.0827 

45 

45.6283 

0.8166 

1.8147% 

1.8147 


Table 50. Results for the OH-6A helieopter with baseline data (HL); network eonfiguration 

14-25-25-1; Ext. DBD learning rule. 
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Figure 55. Results for the OFI-6A helicopter with baseline data (FIL); network configuration 

14-25-25-1; Ext. DBD learning rule, 
c. 14-25-25-1 BPNNExt. DBD (Pruned) 

The Ext. DBD pruned network results are presented because they were 
better than those of the NCD pruned network. The RMS error for this configuration is 
found to be 0.0549. The absolute maximum speed error is 3.43 knots (at 5 knots). The 
maximum airspeed error SD at 1 g is less than 1 knot at all speeds. The overall error SD 
is 0.6505 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.6505 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.3968 

0.3774 

7.5484% 

3.4378 

15 

15.6340 

0.5457 

3.6383% 

1.4552 

25 

23.7417 

0.8713 

3.4852% 

3.3192 

35 

36.3556 

0.6385 

1.8243% 

2.4002 

45 

45.9307 

0.7895 

1.7544% 

2.1794 


Table 51. Results for the OH-6A helicopter with baseline data (HE); network configuration 

14-25-25-1; Ext. DBD (Pruned) learning rule. 
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Figure 56. Results for the OFI-6A helicopter with baseline data (FIL); network configuration 

14-25-25-1; Ext. DBD (Pruned) learning rule. 


7, OH-6A HL Baseline Data IGE Analysis 
a. 14-25-25-1 NCD IGE 

The findings for this configuration are presented in Table 57 and Figure 
57. Results show that the RMS error is 0.0446 and the absolute maximum speed error 
is2.49 knots (at5 knots). The airspeed error SD at 1 g is less than 1 knot at all speeds. The 
overall error SD is 0.6296 knots. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.6296 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.2718 

0.4172 

8.8436% 

2.4936 

15 

15.7052 

0.7943 

5.2950% 

1.9305 

25 

23.5417 

0.5634 

2.2538% 

2.2911 

35 

35.6807 

0.7965 

2.2758% 

1.8238 

45 

45.5805 

0.5918 

1.3150% 

1.4307 


Table 52. Results for the OH-6A helicopter at 12 ft with baseline data (HE); network 

configuration 14-25-25-1; NCD learning rule. 

83 































Predicted Speed vs Grossweight 


Predicted Speed vs Actual Speed 


45 

§35 

"O 

4> 

41 

-25 

■o 

41 

r 

I 15 
□. 

5 

0 


30 

£25 


@@@@@@@@@@@ 

^45 

@ 


2 

35 

• 

eeooooooooo 

& 

i25 

0 


ts 

£ 15 

a. 

• 

@@000000000 

5 

0 

_._._^_ : 


15 

5 

0 


1.5 1.8 2 2.2 

Gross Weight (Iti’^IOOO) 

Slow Speed vs Angle 


2.5 


0 5 


45 


15 25 35 

A^ual Speed (hi) 

Fast Speed vs Angle 


000000®®0000 

< 0 < < < <1 < < ^ < < < 

OOOOOOOOOOOO 


0 60 120 180 240 300 360 

Side Slip Aigle (Deg) 



Side Slip Aigle (Deg) 


14-25-25-1 BPNN NCD QH-6A RMS = 0.0446 ALT= 12 ft tGB IHiah Level! 


Figure 57. Results for the OH-6 A helieopter at 12 ft with baseline data (HL); network 

eonfiguration 14-25-25-1; NCD learning rule. 

b. 14-25-25-1 Ext DBDIGE 

Table 53 and Figure 58 present the findings for this eonfiguration. Results 
show that the RMS error is 0.0682, and the absolute maximum speed error is 4.55 knots, 
whieh oeeurs at 5 knots speed. The maximum airspeed error SD at 1 a is 1.47 knots (at 
25 knots). The overall error SD is 0.9027 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 0.9027 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.8334 

0.4224 

8.4479% 

1.8070 

15 

18.5289 

0.4925 

3.2834% 

4.5586 

25 

24.4931 

1.7556 

5.9025% 

3.3707 

35 

35.5342 

0.8230 

2.3513% 

1.7742 

45 

45.4751 

0.7995 

1.7767% 

1.8916 


Table 53. Results for the OH-6 A helieopter at 12 ft with baseline data (HL); network 

eonfiguration 14-25-25-1; Ext. DBD learning rule. 
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Figure 58. Results for the OH-6 A helicopter at 12 ft with baseline data (HL); network 
configuration 14-25-25-1; Ext. DBD learning rule 
c. 14-25-25-1 NCD (Pruned) IGE 

Figure 59 and Table 54 present the findings for this configuration. Results 
show that the RMS error is 0.0732. The absolute maximum speed error equal to 2.09 
knots occurs at 25 knots speed. The maximum airspeed error SD at 1 g is less than 1 knot 
at all speeds. The overall error SD is 0.6655 knots. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.6655 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.6158 

0.3581 

7.1625% 

2.0768 

15 

15.4433 

0.8653 

5.7688% 

1.9782 

25 

24.1375 

0.6593 

2.6372% 

2.0969 

35 

35.1451 

0.6882 

1.9693% 

1.3938 

45 

44.5488 

0.6975 

1.5501% 

1.4001 


Table 54. Results for the OH-6A helicopter at 12 ft with baseline data (HE); network 

configuration 14-25-25-1; NCD learning rule. 
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Figure 59. Results for the OFI-6A helicopter at 12 ft with baseline data (FIL); network 
configuration 14-25-25-1; NCD (Pruned) learning rule 

8, OH-6A HL Baseline Data OGE Analysis 
a. 14-25-25-1 NCD OGE 

Results for this setup are shown in Figure 60 and Table 55. They show 
that the RMS error is 0.0561. The absolute maximum error equal to 3.15 knots occurs at 
25 knots. The maximum airspeed error SD at 1 g is 1.07 knots, which is observed for 35 


knots speed. The overall error SD is 0.7393 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESUETS 

Total SD = 0.7393 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

3.0288 

0.4230 

8.4591% 

2.7569 

15 

16.1439 

0.7843 

5.2287% 

2.1210 

25 

23.0594 

0.8394 

3.3576% 

3.1582 

35 

35.8068 

1.0779 

3.0798% 

2.5005 

45 

46.0707 

0.5608 

1.2461% 

2.1869 


Table 55. Results for the OH-6A helicopter at 100 ft with baseline data (HL); network 

configuration 14-25-25-1; NCD learning rule. 
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Figure 60. Results for the OH-6A helicopter at 100 ft with baseline data (HL); network 

configuration 14-25-25-1; NCD learning rule. 

b. 14-25-25-1 Ext DBD OGE 

Results are shown in Figure 61 and Table 56. They show that the RMS 
error is 0.0853 and the absolute maximum speed error is 5.39 knots (at 15 knots). The 
maximum airspeed error SD at 1 a equal to 1.67 knots is observed for 25 knots speed. 
The overall error SD is 1.0233 knots. 


Actual 

Airspeed (kts) 

NEURAL NETWORK RESULTS 

Total SD = 1.0233 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 a (kts) 

Percent Error 

at 1 a 

Abs. Maximum 
Error (kts) 

5 

3.7411 

0.4534 

9.0672% 

2.0195 

15 

19.6126 

0.5194 

3.4625% 

5.3927 

25 

25.1916 

1.6786 

6.7143% 

3.7711 

35 

35.3858 

1.070 

3.0573% 

2.0827 

45 

45.7815 

0.8118 

1.8040% 

1.9277 


Table 56. Results for the OH-6A helicopter at 100 ft with baseline data (HL); network 

configuration 14-25-25-1; Ext. DBD learning rule. 
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Figure 61. Results for the OH-6A helicopter at 100 ft with baseline data (FIL); network 
configuration 14-25-25-1; Ext. DBD learning rule, 
c. 14-25-25-1 Ext. DBD (Pruned) OGE 

Figure 62 and Table 57 present the findings for this configuration. Results 
show that the RMS error is 0.0565, which is slightly greater than that of the NCD setup. 
The absolute maximum speed error equal to 4.71 knots occurs at 35 knots speed. The 
maximum airspeed error SD at 1 g is 1.25 knots (at 35 knots). The overall error SD is 
0.8344 knots. 


Actual 

Airspeed (kts) 

NEURAL NI 
Tota 

iTWORK RESULTS 

SD = 0.8344 

Mean of 
Airspeed (kts) 

Airspeed Error 
at 1 G (kts) 

Percent Error 

at 1 G 

Abs. Maximum 
Error (kts) 

5 

2.9652 

0.3962 

7.9235% 

2.7769 

15 

16.1362 

0.6021 

4.0141% 

2.1404 

25 

24.9846 

1.1310 

4.5241% 

2.8046 

35 

37.5287 

1.2565 

3.5785% 

4.7113 

45 

46.6852 

0.7609 

1.6906% 

3.0221 


Table 57. Results for the OH-6A helicopter at 100 ft with baseline data (HE); network 


configuration 14-25-25-1; Ext. DBD (Pruned) learning rule. 
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Figure 62. Results for the OH-6A helicopter at 100 ft with baseline data (FIL); network 
configuration 14-25-25-1; Ext. DBD (Pruned) learning rule. 

C. NETWORK PERFORMANCES SUMMARY 

Results for all experiments are summarized in Table 58 to Table 60. Analysis 
results for the UH-60A helicopter and OH-6A helicopter at sea level and high level 
analyses results are shown separately in Tables 58, 59 and 60 in terms of airspeed error 
SD, absolute maximum error, RMS error and maximum percent error of each NN 
architecture. 


Table 58 displays results obtained for the UH-60A helicopter model. For the UH- 
60A helicopter at 85 ft (out of ground effect condition) a 2-layer BPNN network with the 
NCD learning rule with pruning yielded the best results. This configuration shows a 
predicted airspeed with 0.7 knots error SD, while the maximum error is 3.34 knots, the 
maximum error percentage error is 4.22 % and the RMS error is 0.0355. When the 
helicopter is at 20 ft. (in ground effect condition), the network with NCD rule produced 
an estimate with 0.84 knots error SD, a RMS error equal to 0.0374, a 3 knots maximum 
error and a maximum error percentage within 5.2 %. These results are significantly better 

than those obtained by Haas and McCool, with real flight data. Their study showed that 
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using real flight data as input to the NN, UH-60A airspeed ean be predieted with an 
aeeuracy of ± 5 knots when the aireraft is in ground effeet. However they performed that 
analysis with a referenee airspeed uncertainty of ± 2 knots. Note that prediction accuracy 
improved considerably with the simulator data, which has no wind effect nor any other 


uncertainties. 


UH-60A Helicopter 

BPNN Models 

RMS 

Error 

Error of 

SD at I cj 

Abs. Max. Error 
(knots) 

Max. Percent 

Error (%) 

OGE NCD 

0.0658 

0.7506 

3.85 

(at 35 knots) 

6.4 

(at 5 knots) 

OGE Ext. DBD 

0.0593 

0.7017 

3.86 

(at 45 knots) 

8.49 

(at 5 knots) 

OGE NCD Prune 

0.0355 

0.7167 

3.34 

(at 35 knots) 

4.22 

(at 5 knots) 

Simplified Data OGE NCD 

0.0656 

1.0032 

3.89 

(at 35 knots) 

5.5 

(at 5 knots) 

Simplified Data OGE Ext. DBD 

0.0565 

1.1692 

5.4 

(at 45 knots) 

11.75 

(at 5 knots) 

Simplified Data NCD Prune 

0.0648 

1.1942 

5.4 

(at 25 knots) 

7.95 

(at 5 knots 

One Layer NCD 

0.0407 

0.8394 

3.37 

(at 35 knots) 

4.5 

(at 15 knots) 

One Layer Ext. DBD 

0.0657 

0.7002 

3.79 

(at 35 knots) 

8.34 

(at 5 knots) 

IGE NCD 

0.0374 

0.8469 

3.0 

(at 35 knots) 

5.2 

(at 35 knots) 

IGE Ext. DBD 

0.0637 

0.8308 

4.8 

(at 45 knots) 

6.0 

(at 5 knots) 

IGE NCD Prune 

0.0724 

0.8864 

4.86 

(15 knots) 

5.22 

(at 25 knots) 

Baseline Data NCD 

0.0762 

1.5664 

6.24 

(at 45 knots) 

8.0 

(at 5 knots) 

Baseline Data Ext. DBD 

0.0501 

0.7320 

4.5 

(at 35 knots) 

10.7 

(at 5 knots) 

Baseline Data NCD Prune 

0.1213 

2.0830 

8.64 

(at 45 knots) 

11.5 

(at 25 knots) 

Baseline Data IGE NCD 

0.0759 

1.5516 

6.24 

(at 35 knots) 

10.12 

(at 5 knots) 

Baseline Data IGE Ext. DBD 

0.0499 

0.6505 

4.5 

(at 35 knots) 

8.7 

(at 5 knots) 

Baseline Data IGE NCD Prune 

0.0901 

1.7173 

7.91 

(at 45 knots) 

9.12 

(at 25 knots) 

Baseline Data OGE NCD 

0.0766 

1.5516 

5.86 

(at 25 knots) 

9.3 

(at 25 knots) 

Baseline Data OGE Ext. DBD 

0.0482 

0.6850 

3.8 

(at 35 knots) 

9.5 

(at 5 knots) 

Baseline Data OGE NCD Prune 

0.1087 

1.8642 

7.14 

(at 35 knots) 

9.42 

(at 25 knots) 


Table 58. Overall Results for the UH-60A Helicopter. 
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Results showed that the network with the Ext. DBD rule produeed the best results 
for the baseline data analysis. This network predieted airspeed with a 0.73 knots error SD 
and a 0.05 RMS error. The maximum error was estimated at 4.5 knots at 35 knots speed, 
while the maximum error pereentage was 10.7 % observed for a 5 knots speed. Using the 
baseline data for the IGE eondition airspeed predietion was good with the Ext. DBD 
network sinee it produeed an error SD of 0.65 knots, a RMS error of 0.05, maximum 
error of 4.5 knots at the speed of 35 knots and a maximum pereent error of 8.7%. The 
baseline data OGE eondition results were very elose to those of the baseline data IGE 
eondition results. Using this data, the network with Ext. DBD rule predieted the airspeed 
with a 0.048 RMS error, a 0.68 knots error SD, a 3.8 knots maximum error at the speed 
of 35 knots and a 9.5% maximum pereent error. The single eondition data predietion was 
slightly better than that observed with the baseline data. 

We note that the maximum error and maximum error SD oeeurred mostly at 35 
knots speed for all networks. A potential explanation might be that the hover to forward 
flight translational lift was set to 30 knots for all the EEIGHTLAB simulator models. The 
network was trained with the test data set and tested with the training data set in order to 
explore this idea. Results obtained with the switehed data showed that the maximum error 
and maximum error SD oeeurs when the helieopter is moving at 30 knots. While not 
eonelusive, this set indieates the diffieulty is assoeiated with the simulated helieopter 
performanee near these speeds. 

OH-6A analyses were eondueted at sea level and at high altitude (at 6000 feet 
pressure altitude) using a similar methodology for the UH-60A model. Table 59 shows 
the results of the OH-6A analyses at sea level, and Table 60 shows the results for high 
altitude analyses. 

At sea level, the OGE eondition using the single eondition data NCD network 
with the pruning funetion predieted an airspeed with a 0.76 knots error SD and a 0.067 
RMS error. The maximum error was about 4.69 knots for a 35 knots speed while the 
maximum error was observed at 9.9 %. At sea level, the IGE eondition NCD network 
produeed a 0.84 knots error SD. The RMS error was about 0.059, the maximum error was 
5.1 knots, and the maximum pereent error was 9.7%. Results showed that using baseline 
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data reduced the network performance, as the error SD was 1.169 knots, the RMS error 
was 0.0611, the maximum error was 5 knots and the maximum percent error was 6.8 %. 
The NCD network with the pruning function yielded the best results for this condition. 

For baseline data, the IGE condition NCD pruned network predicted the airspeed 
with a 1.15 knots error SD, while the RMS error was about 0.059 and the maximum error 
was 4.5 for a 35 knots speed. The maximum error percentage was found to be 8.19% for 
this condition. For the baseline data OGE condition, airspeed was estimated with a 1.20 
knots error SD and a 0.065 RMS error with the pruned NCD network. A 5.1 knots 
maximum error occurred for a 35 knots speed and the maximum error percentage was 


6.9%. 


OH-6A Helicopter 

BPNN models (at sea level) 

RMS 

Error 

Error of 

SD at Id 

Abs. Max. 

Error (knots) 

Max. Percent 
Error (%) 

OGE NCD 

0.06 

0.79 

5.2(at 35 knots) 

9.5 (at 5knots) 

OGE Ext. DBD 

0.0735 

1.32 

5 (at 35 knots) 

9.4 (at 5knots) 

OGE NCD Prune 

0.067 

0.759 

4.69(at 35 knots) 

9.9 (at 5knots) 

IGE NCD 

0.059 

0.8465 

5.1 (at 35 knots) 

9.7 (at 5knots) 

IGE Ext. DBD 

0.062 

0.9678 

5.9(at 35 knots) 

12.2 (at 5knots) 

IGE NCD Prune 

0.0674 

0.88 

4.41(at 35 knots) 

12.1 (at 5knots) 

Baseline Data NCD 

0.062 

1.192 

4.82(at 35 knots) 

8.69(at 5knots) 

Baseline Data Ext. DBD 

0.0687 

1.407 

7.44(at 35 knots) 

14(at 5knots) 

Baseline Data NCD prune 

0.0611 

1.16 

5 (at 35 knots) 

6.8(at 5knots) 

Baseline Data IGE NCD 

0.0593 

1.15 

4.49(at 35 knots) 

9.29(at 5knots) 

Baseline Data IGE Ext. DBD 

0.0686 

1.32 

7.44(at 35 knots) 

11.9(at 5knots) 

Baseline Data IGE NCD Prune 

0.0592 

1.15 

4.51 (at 35 knots) 

8.19(at 5knots) 

Baseline Data OGE NCD 

0.0671 

1.20 

4.82(at 35 knots) 

7.2(at 5knots) 

Baseline Data OGE Ext. DBD 

0.0687 

1.3 

6.6(at 35 knots) 

9.4(at 5knots) 

Baseline Data OGE NCD Prune 

0.0650 

1.2 

5.1 (at 35 knots) 

6.9(at 5knots) 


Table 59. Overall Results for OH-6A Helicopter at Sea Eevel. 

The pruned NCD network predicted the airspeed with a 0.64 knots error SD for 
the OH-6A high level OGE condition. The RMS error was found to be equal to 0.0475, 
while the maximum error was 2.8 knots at 5 knots speed and the maximum error was 
about 7%. Einally, the airspeed was predicted with a 0.64 knots error SD, a 0.0433 RMS 
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error, a 2.58 knots maximum error and a 5.6% maximum error pereentage for the IGE 
eondition using the single eondition data with the pruned NCD network. 

The following results were obtained for the baseline data. The Ext. DBD network 
performed best with a 0.65 knots error SD, a 0.055 RMS error, a 3.4 knots maximum 
error for a 5 knots speed and a 7.5 % maximum error pereentage. The NCD pruned 
network yielded the best results for the baseline data IGE eondition. Eor this setup, the 
RMS error was 0.0352 while the error SD was 0.66 and the maximum error was 2 knots 
at the speed of 25 knots. The maximum pereent error was 7.1% for this eondition. The 
baseline data OGE eondition results were elose to the IGE eondition results. The NCD 
network estimated the airspeed with an error SD of 0.74 knots, a RMS error of 0.056, a 
maximum error of 3.15 knots and a maximum pereent error of 8.4%. 


OH-6A Helicopter 

RMS 

Error of 

Ahs. Max. 

Max. Percent 

BPNN models (at high altitude) 

Error 

SD at la 

Error (knots) 

Error (%) 

OGE NCD 

0.0485 

0.6637 

2.47 (at 5knots) 

9.47(at 5knots) 

OGE Ext. DBD 

0.0958 

1.022 

6.36(at 15knots) 

5.67(at 5knots) 

OGE NCD Prune 

0.0475 

0.64 

2.85(at 5knots) 

7. Rat 5knots) 

IGE NCD 

0.0516 

1.20 

2.59(at 5knots) 

9.12(at 5knots) 

IGE Ext. DBD 

0.0865 

0.95 

5.8(at 15knots) 

5.53(at 5knots) 

IGE NCD Prune 

0.0443 

0.64 

2.58(at 5knots) 

5.59(at 5knots) 

Baseline Data NCD 

0.05 

0.71 

3.15(at 25knots) 

8.73(at 5knots) 

Baseline Data Ext. DBD 

0.077 

1.02 

5.4(at 15knots) 

8.8(at 5knots) 

Baseline Data Ext. DBD prune 

0.055 

0.65 

3.43(at 5knots) 

7.5(at 5knots) 

Baseline Data IGE NCD 

0.0446 

0.63 

2.49(at 5knots) 

8.84(at 5knots) 

Baseline Data IGE Ext. DBD 

0.0687 

0.9 

4.55(at 5knots) 

8.44(at 5knots) 

Baseline Data IGE NCD Prune 

0.0352 

0.66 

2.09(at 25knots) 

7. Rat 5knots) 

Baseline Data OGE NCD 

0.056 

0.74 

3.15(at 25knots) 

8.45(at 5knots) 

Baseline Data OGE Ext. DBD 

0.0853 

1.02 

5.39(at 15knots) 

9(at 5knots) 

Baseline Data OGE Ext. DBD Prune 

0.0565 

0.83 

4.7 Rat 35knots) 

7.9(at 5knots) 


Table 60. Overall Results of OH-6A Helieopter at High Altitude. 

In summary, these results show that the NN approaeh to prediet airspeed using 
simulation data is quite promising. The BPNN with two hidden layers and 25 PEs in eaeh 
layer performs the best among all studied arehiteetures. We note that different learning 
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rules yielded different results and that enabling the pruning faeility improved the network 
performanee in most oases. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

Today, military helicopters perform a wide variety of tasks in conditions ranging 
from hot and dry to cold and wet, windy and low visibility weather. Accurate low speed 
velocity sensing devices are essential because aircraft velocity and position information 
are what pilots need to perform safely in these regimes. However, conventional speed 
measuring systems do not work accurately when the aircraft speed is below 40 knots. 

NN-based airspeed prediction studies developed by McCool, Haas and others 
showed that NN is a good approach to solve this problem. The objective of this thesis is 
to build a NN model in order to predict OH-6A helicopter airspeed in the low speed 
environment using a flight simulator to obtain the parameters required for the NN. In our 
study the NN model is developed with Neural Works Professional Plus/II software. 

First a NN model of the UH-60A helicopter was built and implemented, and 
several NN configurations analyzed. The reason for building a NN model for the UH- 
60A helicopter was to lay out a background to make a comparison of NN predictions 
using simulator data and NN predictions using real flight data. The results showed using 
simulator data potentially improves the accuracy of prediction significantly. 

Three different methods were investigated to select the NN training data. The first 
one is a single condition data set in which the data belongs to one altitude only. The 
second one is called the baseline data set, and is formed by combining the data of two 
single condition data sets of different altitudes. The third set is obtained by applying 
principal component analyses to decrease the input space dimension and it is called a 
simplified data set. Results showed that the network trained by using a single condition 
data set proved to be the most successful and performance degraded only slightly with the 
baseline data. Moreover, the BPNN network produced more successful predictions than 
the RBFN implementation. 

Among all BPNN architecture types considered, a two-hidden layer BPNN with 
an enabled pruning facility for the NCD learning rule showed the best performance. At 
sea level pressure altitude, the UH-60A low airspeed was predicted with one-sigma 
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accuracy of ±0.71 knots when the aireraft was out of ground effeet. The aeeuraey of 
predietion was ±0.88 knots when the aireraft was in ground effeet. 

The OH-6A low speed was predieted using a similar methodology based on the 
results obtained from the UH-60A model. The OH-6A analysis was performed at high 
pressure altitude as well as at seal level altitude. Results showed that at sea level, the OH- 
6A airspeed eould be predieted with one-sigma aeeuraey of ±0.75 knots when the 
aireraft is out of ground effeet. The best performanee for all high altitude analyses 
obtained from the pruned NCD network was with single eondition data. For IGE 
conditions, the predietion aeeuraey was about ±0.88 knots. High altitude analysis was 
performed at 6000 feet. The results showed that at high altitude, the OH-6A airspeed 
eould be predieted with an aeeuraey of ± 0.64 knots when the aireraft is out of ground 
effeet. For IGF eonditions, the predietion aeeuraey was about ± 0.64 knots. 

This study showed that a NN based approaeh to determine OH-6A helieopter 
airspeeds using a flight simulator is quite promising. The approaeh considered presents a 
meehanieally simple alternative to eurrent low airspeed measurement systems, and as a 
result eontributes to inerease the flight safety and eombat effeetiveness. 

B, RECOMMENDATIONS FOR FURTHER RESEARCH 

Several avenues are available to extend the work presented in this study. First, 
results obtained using real flight data are needed. These results should be eompared to 
the results obtained here to evaluate the performanee of our approaeh. Seeond, more 
sophistieated maneuvers need to be simulated and analyzed. The NN model might be 
further improved by using different network models, model parameters, ete. Another type 
of network to be investigated may be NN-based sideslip angle estimator, whieh would 
provide aeeurate wind direetion for pilot-aid systems. Finally, the effeets of measurement 
noise, sensor errors, and sensor failure should be investigated. 
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APPENDIX A. NEURAL WORKS PROFESSIONAL PLUS/II 

PROGRAM SETUP 


Neural Works Professional II/PLUS is a powerful and flexible development tool. 
It has over 20 different algorithms ineluding eornmon network types sueh as Baek 
Propagation, Kohonen and Radieal Basis Funetions. The produet ineludes a variety of 
diagnostie tools as well as options for quiek and easy building networks. 

Before ereating the network, data files must be put under the direetory of 
Neuralware Professional direetory. The test and training input data files must be given 
different names but both of them must have extension of “.nna”. The data files used in 
this thesis may be obtained from Professor R. W. Duren, Naval Post Graduate Sehool 
Department of Aerospaee and Aeronautieal Engineering, by request. 

Neural network development proeess begins with eolleeting and preproeessing 
input data. In this study data is eolleeted from FLIGHTLAB simulation and prepared by 
using MATLAB. Preparation involves eneoding data to a format that NN ean deal with. 
For this problem data set was eonverted to a matrix with a dimension of 1114x14, where 
column s represent the inputs and rows represent the samples of each inputs. MATFAB is 
also used to separate the data into to sets, training set and test set, as well as to display the 
results. MATFAB codes for preprocessing and post processing the data are included in 
Appendix B. 

After starting NeuralWorks program. Back-propagation command, under the 
InstaNet menu, must be selected to create a BP network. Selecting this command opens 
up BP dialog box, which allows user to build the main frame of the network by entering 
number of layers and number of PEs per layer. In addition to these, learning coefficient 
of each layer, momentum term, learning rule, transfer function, test data and training data 
of the network can be selected. Also by selecting “minmax” radio button, inputs can be 
mapped from the data file to a desired range, such as -1 to 1. 

Activating this box opens up another window, which asks user how to display the 
performance of the network, such as RMS error, classification rate, etc. After these 
selections network is created which is presented in a window like in Figure 5. Created 
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network can be saved either in binary format or in ASCII format by selecting the Save 
command under the File menu. 

Selecting the Learn command under the Run menu starts training. At this point, 
number of learning iterations can be entered by the user. As each training example is 
presented to the network, the network produces an output, which is used to evaluate the 
training performance of the network. Another way to train the network is to use Savebest 
command under the Run menu. This commands opens up Run/Check dialog box, which 
makes pruning facility accessible. Based on a decision criteria specified by the user, 
pruning facility disables connections in a network as the network is training. 

Network is tested using the test data sets and by selecting the Test command 
under the Run menu. The desired outputs, along with the actual network results are 
written to the results file which has “.nnr” extension. 


Based on the above explanations and after preparing the data, OH-6A helicopter 
14-25-25-1 BPNN NCD (Prune) network model created using the following steps: 

Start Neural Works on the computer 

Select Back-Propagation command from the InstaNet menu. This command pops 


up the following window. 


NeuralWorks Professional II/PLU5 


File InstaNet I/O Instrument Run 
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InstaNet / Back Propagation 
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Set Epoch From File 


ExtDBD 
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Delta-Bar-Delta 


Learn 
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Sine 


Rcl/Test 


xj 


|quad_tst 

|quad_tst 

lauad trn.nnc.^1 

Iquad trn.nnc.^1 

recirc.nna 

recurtrn.nna_1 

recurtst.nna 
sample, nna ^ 

recirc.nna 
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OK I Cancel | 


Help 


NeuralWorKs Professional _ . , , 

Copvrighi2000 NeuralWare. All rights reserved. 


Figure 63. Back-propagation network setup window 
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1. Select “trainfile tm” from the training input scroll window. 

2. Select “testfile tst” from the testing input scroll window. 

3. In the number of PEs section, enter the following numbers: 

-Input: 14 
-Hidl:25 
-Hid2: 25 
-Output: I 

4. Enter 0.4 for momentum and 0.5 for LCoef ratio. 

5. Select Norm-Cum-Delta for the learning rule. 

6. Select TanH for the transfer function. 

7. Check MinMax table box. 

8. Click the OK button. 

9. After clicking OK the following window opens automatically. Select RMS 
Error. 



Eigure 64. Instrument /Create menu 

10. Click OK. 

11. Start training by selecting the SaveBest command under the Run menu. 
SaveBest command setup is shown in Eigure #. Training can also be started by 
selecting Learn command under the same menu. SaveBest command allows 
user to use Pruning facility. 

12. Type the name of the file that results are to be written. 

13. Enter 50000 in the Eor field. 

14. Enter 1000 for the Test Interval. 
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15. Check the pruning radio button. 

16. Enter 0.975 for the Tolerance field. 

17. Select Classification Rate for the Objective Function list. 

18. Click OK. to start training. 



Figure 65. SaveBest command window 
19. After training is completed select Test in the Run menu. 



Figure 66. Test command window 

20. Select One Pass/All. 

21. Click OK. 

After test process has been completed the performance of the network is displayed 
on the main window in terms of selected options, such as RMS error, network weights, 
etc. and the results are stored in the test file with .nnr extension. In this work, the result 
files were saved as text files in order to be exported into the MATFAB. 
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APPENDIX B. MATLAB^ M-FILES 


MATLAB m-files were developed to prepare the simulation outputs for NN and 
to process the NN results. 


a. FILE Name: Gregcode 

The first code included in this appendix is the FLIGHLAB scope language, an 
interpretive language that uses the industry standard MATLAB syntax. This routine 
loads the specified model and runs the simulator. The outputs of this routine are used as 
inputs to the NN model. 

This code was developed by LT. Gregory OUELLETTE, USN. for the Naval 
Postgraduate School in partial fulfillment of the requirements for a Masters Degree in 
Aeronautical Engineering. 

September 2001. 

exec("$EE_DIR/flme/models/articulated/arti-rgd-3iv-qs-simeng.def',l) 
world_model_airframe_cpg_testcond_poszic = -85; 
exec("xatestcond.exc",l) 
exec("xamodeltrim.exc",l) 

goto world 
group test 

gw = [16000:1000:24000]’; 
hdg = [30:30:360]'; 
vel= [0:5:30]’; 
column = 0; 

utrim = @trimvariable; 

statesave = savestates(world_topsolve); 

for ngw = 1 :prod(size(gw)) 
for nhdg = 1: prod(size(hdg)) 
for nvel = 1: prod(size(vel)) 

world_model_airframe_cpg_testcond_poszic = -85; 
worldmodelairframecpgtestcondveq = vel(nvel); 
world_model_airframe_cpg_testcond_gamh = hdg(nhdg); 
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world_model_data_vweight = gw(ngw); 

exec("xatestcond.exc",l) 

exec("xamodeltrim.exc",l) 

outputs(column+nvel,l) = world_model_control_data_xatrm; 
outputs(column+nvel,2) = world_model_control_data_xbtrm; 
outputs(column+nvel,3) = world_model_control_data_xctrm; 
outputs(column+nvel,4) = world_model_control_data_xptrm; 
outputs(column+nvel,5) = world_model_airframe_cpg_xaout_p; 
outputs(column+nvel,6) = world_model_airframe_cpg_xaout_q; 
outputs(column+nvel,7) = world_model_airframe_cpg_xaout_r; 
outputs(column+nvel,8) = world_model_airframe_cpg_xaout_phi; 
outputs(column+nvel,9) = world_model_airframe_cpg_xaout_psi; 
outputs(column+nvel, 10) = world_model_airframe_cpg_xaout_radralt; 
outputs(column+nvel,l 1) = world_model_airframe_cpg_xaout_vclimb; 
outputs(column+nvel, 12) = world_model_rotor l_rotor_cpg_xaout_omega; 
outputs(column+nvel,13) = world_model_propulsion_cpg_xaout_etorq; 
outputs(column+nvel,14) = world_model_data_vweight; 
outputs(column+nvel,15) = world_model_airframe_cpg_xaout_tas; 
outputs(column+nvel, 16) = world_model_airframe_cpg_testcond_gamh; 
outputs(column+nvel,17) = world_model_airframe_cpg_testcond_veq; 
outputs(column+nvel, 18) = 0; 

end 

column = column + prod(size(vel)); 
savestates(world_topsolve,statesave); 

@trimvariable = utrim; 

end 

end 

trainset = outputs'; 
save("velhead.sav",tramset) 


exec("$FL_DlR/flme/models/articulated/arti-rgd-3iv-qs-simeng.def',l) 
world_model_airframe_cpg_testcond_poszic = -85; 
exec("xatestcond.exc",l) 
exec("xamodeltrim.exc",l) 

goto world 
group test 
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gw = [16000:1000:24000]’; 
hdg = [-60:30:60]’; 
vel= [35:5:40]’; 
column = 0; 

utrim = @trimvariable; 

statesave = savestates(world_topsolve); 

for ngw = 1 :prod(size(gw)) 
for nhdg = 1: prod(size(hdg)) 
for nvel = 1: prod(size(vel)) 

world_model_airframe_cpg_testcond_poszic = -85; 
worldmodelairframecpgtestcondveq = vel(nvel); 
world_model_airframe_cpg_testcond_gamli = hdg(nhdg); 
world_model_data_vweight = gw(ngw); 

exec("xatestcond.exc",l) 

exec("xamodeltrim.exc",l) 

outputs(column-l-nvel,l) = worldmodelcontroldataxatrm; 
outputs(column+nvel,2) = worldmodelcontroldataxbtrm; 
outputs(column-l-nvel,3) = worldmodelcontroldataxctrm; 
outputs(column-l-nvel,4) = worldmodelcontroldataxptrm; 
outputs(column-l-nvel,5) = world_model_airframe_cpg_xaout_p; 
outputs(column+nvel,6) = world_model_airframe_cpg_xaout_q; 
outputs(column-l-nvel,7) = world_model_airframe_cpg_xaout_r; 
outputs(column-l-nvel,8) = world_model_airframe_cpg_xaout_phi; 
outputs(column-l-nvel,9) = world_model_airframe_cpg_xaout_psi; 
outputs(column+nvel, 10) = world_model_airframe_cpg_xaout_radralt; 
outputs(column-l-nvel,l 1) = world_model_airframe_cpg_xaout_vclimb; 
outputs(column-l-nvel, 12) = world_model_rotor l_rotor_cpg_xaout_omega; 
outputs(column+nvel,13) = world_model_propulsion_cpg_xaout_etorq; 
outputs(column-l-nvel,14) = world_model_data_vweight; 
outputs(column-l-nvel,15) = world_model_airframe_cpg_xaout_tas; 
outputs(column-l-nvel, 16) = world_model_airframe_cpg_testcond_gamh; 
outputs(column+nvel,17) = world_model_airframe_cpg_testcond_veq; 
outputs(column+nvel, 18) = 0; 
end 

column = column + prod(size(vel)); 
savestates(world_topsolve,statesave); 

@trimvariable = utrim; 

end 

end 
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trainset = outputs'; 

save(" velheadfast. sav",trainset) 


b, FILE Name: Ozcan.m 

% This routine takes simulation output files as input and makes them applieable for NN 
% model by forming a matrix, where eolumns represent the input parameters to NN and 
% rows represent the samples. 

% matl is the output matrix of the routine where all the data are stored and to be used as 
% input data to the NN. 

% The last part of the routine may be used to obtain the simplified data using 
% eigenvalues and eigenveetors. Dataeig is the resultant matrix for simplified data. 

% September 2001 

elear all; 
format short e 
load velhead.txt; 
load velheadfast.txt; 

m=l; 

for i = l;length(velhead)/6 

slat_stiok(i) = velhead(m,l); 
slong_stiok(i) = velhead(m,2); 
sooll_pos(i) = velhead(m,3); 
sped_pos(i) = velhead(m+l,l); 
sroll_rate(i) = velhead(m+l,2); 
spitoh_rate(i) = velhead(m+l,3); 
syaw_rate(i) = velhead(m+2,l); 
spitoh_att(i) = velhead(m+2,2); 
sroll_att(i) = velhead(m+2,3); 
salt(i) = velhead(m+3,l); 
solimb_rate(i) = velhead(m+3,2); 
smrb_rpm(i) = velhead(m+3,3); 
seng_torque(i) = velhead(m+4,l); 
sgw(i) = velhead(m+4,2); 

stas_trim(i) = velhead(m+4,3)*(360/608); % ft/see is eonverted into knots 

sheading(i) = velhead(m+5,l); 
stas_trgt(i) = velhead(m+5,2); 


% loading aseii file (simulation output file) 
%loading aseii file 
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szero(i) = velhead(m+5,3); 

m=m+6; 

end 

m=l; 

for i = l;length(velheadfast)/6 

flat_stick(i) = velheadfast(m,l); 
flong_stick(i) = velheadfast(m,2); 
fcoll_pos(i) = velheadfast(m,3); 
fped_pos(i) = velheadfast(m+l,l); 
froll_rate(i) = velheadfast(m+l,2); 
fpitch_rate(i) = velheadfast(m+l,3); 
fyaw_rate(i) = velheadfast(m+2,l); 
fpitch_att(i) = velheadfast(m+2,2); 
froll_att(i) = velheadfast(m+2,3); 
falt(i) = velheadfast(m+3,l); 
fclimb_rate(i) = velheadfast(m+3,2); 
fmrb_rpm(i) = velheadfast(m+3,3); 
feng_torque(i) = velheadfast(m+4,l); 
fgw(i) = velheadfast(m+4,2); 

ftas_trim(i) = velheadfast(m+4,3)*(360/608); % ft/sec is converted into knots 

fheading(i) = velheadfast(m+5,1); 

ftas_trgt(i) = velheadfast(m+5,2); 

fzero(i) = velheadfast(m+5,3); 

m=m+6; 

end 

lat_stick=[slat_stick,flat_stick]; 

long_stick=[slong_stick,flong_stick]; 

coll_pos=[scoll_pos,fcoll_pos]; 

ped_pos=[sped_pos,fped_pos]; 

roll_rate=[sroll_rate,froll_rate]; 

pitch_rate=[spitch_rate,fpitch_rate]; 

yaw_rate=[syaw_rate,fyaw_rate]; 

pitch_att=[spitch_att,fpitch_att]; 

roll_att=[sroll_att,froll_att]; 

alt=[salt,falt]; 

climb_rate= [sclimb_rate,fclimb_rate]; 
mrb_rpm=[smrb_rpm,fmrb_rpm]; 
eng_torque=[seng_torque,feng_torque]; 
gw=[sgw,fgw]; 

tas_trim= [stas_trim,ftas_trim]; 
heading=[sheading,fheading]; 
tas_trgt=[stas_trgt,ftas_trgt]; 

%***for UH-60A helicopter add parameter “eng torque” to the matrix below**** 
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matl=[lat_stick;long_stick;coll_pos;ped_pos;roll_rate;pitch_rate;yaw_rate;pitch_ 

att;roll_att;alt;climb_rate;mrb_rpm;gw;heading;tas_trgt]' 

% To simplify the data by eigenvalues and eigenvectors use the following part 
% mat2=mat 1' * mat 1; 

% [v,d]=eig(mat2); 

% After examining the eigenvalues create submatrix u 
% u=[matl(:,9),matl(:,10),matl(;,ll),matl(:,12),matl(:,13),matl(:,14)]; 

% dataeig=matr*u; 


c, FILE Name: Train.m 

% This routine takes the output matrix of Ozcan.m code, mail or dataeig, by loading the 
% ascii file and filters it so that train data set can be obtained from the whole data. 

% September 2001 

clear all; 
format short e 

load oh6sl_100ft.txt; % loading ascii file 

[rmax,cmax]=size(oh6sl_l 00ft); 

dmax=924 ; % Enter the # of rows associated with slow velocity data 

s=l; 

x=0; 

y=8; 

for d=l:dmax/14 
for i=l;2;7 

trndt(s,:)=oh6sl_100ft(i+x,:); 

s=s+l; 

end 

for j=0:2;6 

trndt(s,:)=oh6sl_100ft(j+y,:); 

s=s+l; 

end 

x=x+14; 

y=y+14; 

end 

for dl=(dmax+2):2;rmax 

trndt(s,:)=oh6sl_100ft(dl,:); 

s=s+l; 

end 

trndt % training data set 
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d. FILE Name: Test.m 


% This routine takes the output matrix of Ozean.m code, mail or dataeig, by loading the 
% ascii file and filters it so that test data set can be obtained from the whole data. 

% September 2001 

clear all; 
format short e 

load oh6sl_100ft.txt; % loading ascii file 

[rmax,cmax]=size(oh6sl_l 00ft); 

dmax=924 ; % Enter the #of rows associated with slow velocity data 

sl=l; 

xl=0; 

yl=9; 

for d2= 1; dmax/14 
for 1=2:2;6 

tstdt(sl,;)=oh6sl_100ft(l+xl,:); 

sl=sl+l; 

end 

for jl=0:2:4 

tstdt(sl,:)=oh6sl_100ft(jl+yl,:); 

sl=sl+l; 

end 

xl=xl+14; 

yl=yl+14; 

end 

for d3=(dmax+l):2:rmax 
tstdt(sl,:)=oh6sl_100ft(d3,:); 
sl=sl+l; 
end 

tstdt % test data set 
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e. FILE Name: Bersan.m 


% This program used to process the outputs of the NN. In this routine the output file of 
% the NN is taken as input and vectors of NN predicted speeds are created related to each 
% gross weight and sideslip angle of the helicopter. This fide also produces the figures 
% and the evaluation of the NN results. 

% October 2001 

% Note; For baseline data and single data use the specified sections of the program. Also 
% depending on the helicopter type activate and deactivate the stated lines of the code. 

clear all; 

load oh6_sl_ncdprune.txt; % loading ascii file 

a=oh6_sl_ncdprune(:,2); % NN outputs(predicted speeds) 

ilk=oh6_sl_ncdprune(:,l); % target values(actual speeds) 

%set initial variables 

g=i; 

j=i; 

k=l; 

1 = 1 ; 

m=l; 

n=length(a); 

aci=[30;30:360]; % angles for slow speed 

acil=[-60;30;60]; % angles for fast speed 

ai=length(aci); 

% x=l.6:0.1:2.4; % gross weight range of uh-60a helicopter(10M lb) 

x=l.500:. 100:2.550; % gross weight range of oh-6a helicopter(10^3 lb) 

%(forUH-60A; 16000; 1000:24000 - for OH-6A; 1500:100:2550 ) 

gw=length(x); 

na=(ai*gw*6); % slow velocity vector dimension (for baseline data 

% multiply by 6, for single data multiply by 3) 

% Classify the speeds and get each speed vectors 
for i=l;3;na 

al(j,l)=a(i,l); %vector of 5 kts 

fark5(j)= abs(al(j,l)-ilk(i,l)); 

j=j+i; 

end 


no 



for i=2;3;na 

a2(k,l)=a(i,l); %vector 15 kts 

fark 15 (k)= abs(a2(k, 1 )-ilk(i, 1)); 
k=k+l; 
end 

for i=3:3;na 

a3(l,l)=a(i,l); %vector25kts 

fark25(l)= abs(a3(l,l)-ilk(i,l)); 

1 = 1 + 1 ; 

end 

for t=(na+l):2;n 

a4(m,l)=a(t,l); % veetor 35 kts 

fark35(m)= abs(a4(m,l)-ilk(t,l)); 
m=m+l; 
end 

for t=(na+2):2;n 

a5(g,l)=a(t,l); %veotor45kts 

fark45(g)= abs(a5(g,l)-ilk(t,l)); 
g=g+i; 
end 

% elassify speeds aceording to gross weights 

y=length(al); 

b=y/gw; 

% for speed = 5 kts 
c=l; 

for i=l;gw 
ifi>=2 
z=l; 

for j=(((i-l)*b)+l):i*b 
all(z,i)=alO‘,l); 
z=z+l; 
end 
else 

for v=l:b 
all(c,l)=al(v,l); 
o=c+l; 
end 
end 
end 

figure(l) 
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subplot(2,2,l);plot(x,al 1,'ro') 

%for speed = 15 kts 
c=l; 

for i=l;gw 
ifi>=2 
z=l; 

for j=(((i-l)*b)+l):i*b 
a22(z,i)=a20‘,l); 
z=z+l; 
end 
else 

for v=l;b 
a22(c,l)=a2(v,l); 
c=c+l; 
end 
end 
end 

hold on 

subplot(2,2,l);plot(x,a22,'g<') 

%for speed = 25 kts 
e=l; 

for i=l;gw 
ifi>=2 
z=l; 

for j=(((i-l)*b)+l):i*b 
a33(z,i)=a3(j,l); 

z=z+l; 

end 

else 

for v=l;b 
a33(c,l)=a3(v,l); 

0 = 0 + 1 ; 

end 

end 

end 

hold on 

subplot(2,2,1 );plot(x,a33 ,'bo') 

% for airspeed = 35 kts 
yl=length(a4); 
bl=yl/gw; 
c=l; 

for i=l:gw 
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ifi>=2 

z=l; 

for j=(((i-l)*bl)+l):i*bl 
a44(z,i)=a40‘,l); 

z=z+l; 

end 

else 

for v=l:bl 
a44(o,l)=a4(v,l); 

0 = 0 + 1 ; 

end 

end 

end 

hold on 

subplot(2,2,l);plot(x,a44,'m>') 

% for airspeed 45 kts 
0 = 1 ; 

for 1=1: gw 
ifi>=2 
z=l; 

for j=(((i-l)*bl)+l):i*bl 
a55(z,i)=a50,l); 

z=z+l; 

end 

else 

for v=l:bl 
a55(o,l)=a5(v,l); 

0 = 0 + 1 ; 

end 

end 

end 

hold on 


%for UH-60A 
% for OH-6A 

%for UH-60A 
% for OH-6A 


% Plotting predioted speed vs aotual speed 
subplot(2,2,2);plot(ilk,a,'o') 


subplot(2,2,1 );plot(x,a55 ,'ko') 
axis([1.4 2.6 0 50]) 

title('Predioted Speed vs Grossweighf,'FontSize',10) 
xlabel('Gross Weight (lb*10000)VFontSize',7) 

% xlabel('Gross Weight (lb*1000)VFontSize',7) 
ylabel('Predioted Speed (Kt)VFontSize',7) 
set(goa,'xtiok’,[1.6 1.8 2.0 2.2 2.4],'ytiok',[0 5 15 25 35 45]) 

% set(goa,'xtiok',[1.5 1.8 2.0 2.2 2.5],'ytiok',[0 5 15 25 35 45]) 
hold off 
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axis([0 55 0 55]) 

title('Predicted Speed vs Aetual SpeedVFontSize',10) 
xlabel('Aetual Speed (Kt)VFontSize',7) 
ylabel('Predieted Speed (Kt)VFontSize',7) 
set(gea,'xtiek',[0 5 15 25 35 45],'ytiek',[0 5 15 25 35 45]) 
hold off 

% plotting airspeed vs. angle (for single data only) 

% subplot(2,2,3);plot(aoi,al 1,'ro') 

% hold on 

% subplot(2,2,3);plot(aoi,a22,'g<') 

% hold on 

% subplot(2,2,3);plot(aoi,a33,'bo') 

% axis([0 400 0 35]) 

% title('Slow Speed vs AngleVFontSize', 10) 

% xlabel('Side Slip Angle (Deg)VFontSize',7) 

% ylabel('Predieted Speed (Kt)VFontSize',7) 

% set(gea,'xtiek',[0 60 120 180 240 300 360],'ytiek',[0 5 15 25 30]) 

% hold off 

% subplot(2,2,4);plot(aoil ,a44,'m>') 

% hold on 

% subplot(2,2,4);plot(aoil,a55,'ko') 

% axis([-100 100 30 50]) 

% title('Fast Speed vs Angle VFontSize', 10) 

% xlabel('Side Slip Angle (Deg)','FontSize',7) 

% ylabel('Predieted Speed (Kt)','FontSize',7) 

% set(gea,'xtiek',[-60 -30 0 30 60],'ytiek',[30 35 40 45 50]) 

% hold off 

% For Baseline data plotting speed vs. angle use the following seetion of the routine 

for df=l;12 

alml(df,:)=al l(df,;); 
a2ml(df,:)=a22(df,:); 
a3ml(df,:)=a33(df,;); 
end 

for vo=l;5 

a4ml(vo,:)=a44(vo,:); 

a5ml(vo,;)=a55(vo,;); 

end 

nd=l; 

for gf= 13:24 

alm2(nd,:)=al l(gf,:); 
a2m2(nd,: )=a22(gf,:); 
a3m2(nd,:)=a33(gf,:); 
nd=nd+l; 
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end 

ndl=l; 

for gfl=6;10 

a4m2(ndl, ;)=a44(gfl, 
a5m2(ndl,;)=a55(gfl,:); 
ndl=ndl+l; 
end 

subplot(2,2,3);plot(aei,alml,'ro') 
hold on 

subplot(2,2,3);plot(aci,alm2,'ro') 
hold on 

subplot(2,2,3);plot(aci,a2ml,'g<') 
hold on 

subplot(2,2,3);plot(aoi,a2m2,'g<') 
hold on 

subplot(2,2,3);plot(aoi,a3ml,'bo') 
hold on 

subplot(2,2,3);plot(aci,a3m2,'bo') 
axis([0 400 0 35]) 

title('Slow Speed vs AngleVFontSize',10) 

xlabel('Side Slip Angle (Deg)VFontSize',7) 

ylabel('Predicted Speed (Kt)VFontSize',7) 

set(gca,’xtiek',[0 60 120 180 240 300 360],'ytiek',[0 5 15 25 30]) 

hold off 

subplot(2,2,4);plot(acil ,a4ml ,'m>') 
hold on 

subplot(2,2,4);plot(acil,a4m2,'m>') 
hold on 

subplot(2,2,4);plot(acil ,a5ml ,'ko') 
hold on 

subplot(2,2,4);plot(aoil,a5m2,'ko') 

axis([-100 100 30 50]) 

title('Fast Speed vs AngleVFontSize', 10) 

xlabel('Side Slip Angle (Deg)VFontSize',7) 

ylabel('Predicted Speed (Kt)VFontSize',7) 

set(gca,'xtiek',[-60 -30 0 30 60],'ytiek',[30 35 40 45 50]) 

hold off 

% Finding max errors for the worst ease 
for i=l:n 

fark(i)=(abs(ilk(i)-a(i))); 

end 

maxfark=max(fark) 
maxfark5=max(fark5 ) 
maxfarkl 5=max(farkl 5) 
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maxfark25=max(fark25) 
maxfarkS 5=max(fark3 5 ) 
maxfark45=max(fark45) 


% Computation of SD and percent error for each speed; 
mn5=mean(al) 
sd5=std(al) 
per5=(sd5/5)*100 
mnl5=mean(a2) 
sdl5=std(a2) 
perl5=(sdl5/15)*100 
mn25=mean(a3) 
sd25=std(a3) 
per25=(sd25/25)*100 
mn35=mean(a4) 
sd35=std(a4) 
per35=(sd35/35)*100 
mn45=mean(a5) 
sd45=std(a5) 
per45=(sd45/45)*100 

% SD for whole data set 

sdw=sqrt((((na/3)-l)*(sd5^2+sdl5^2+sd25^2)+(((n-na)/2)-l)* sd35^2+sd45^2)) /(n-1)) 
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